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ABSTRACT 


DYNAMIC  BEHAVIOR  OF  A  MENISCUS  SUBJECTED  TO  A  TRANSIENT 
ACCELERATION  FIELD 

Name:  Yeikes,  Kirk,  Lee 
University  of  Dayton,  1994 

Advisor:  Dr.  Kevin  P.  Hallinan 

The  dynamic  response  of  an  unheated  and  heated  meniscus  within  a  capiUaiy  tube 
subjected  to  single  and  multiple  cycle  transient  acceleration  induced  forces  was  studied 
experimentally  and  anal)rtically.  A  capillary  tube  was  partially  £Qled  with  either  water 
or  ethyl  alcohol  and  mounted  on  a  centrifuge  for  visual  observation  of  the  dynamic 
response  of  the  meniscus  subject  to  an  acceleration  transient.  A  one-dimensional 
mathematical  model  was  formulated  to  predict  the  djmiamic  response  of  the  meniscus  and 
compared  with  experimental  results.  For  the  case  in  which  the  meniscus  was  heated, 
thermocapiUary  stresses  in  the  near  contact  line  region  were  accoimted  for  in  the  model 
by  a  change  in  contact  angle. 

Comparison  of  the  experimental  results  and  analytical  calculations  for  the  case  of 
an  unheated  meniscus  showed  that  a  combined  axial  and  elevated  transverse  acceleration 
significantly  altered  the  dynamic  response  of  the  meniscus  and  resulting  in  an  increased 
meniscus  velocity  and  subsequently  lower  meniscus  position.  Systems  characterized  by 
large  Bond  numbers,  greater  than  1.5,  as  distinguished  by  the  larger  capillary  tube 

iii 


diameters,  were  particularly  sensitive  to  an  elevated  transverse  acceleration  component 
greater  than  2g.  Systems  typified  by  a  low  Bond  number,  less  than  1.5,  tended  to  be 
insensitive  to  a  transverse  acceleration  component  due  to  the  small  capillary  tube 
diameter.  For  systems  typified  by  the  combination  of  a  low  Bond  number,  less  than  1.5, 
and  high  c^illary  number,  greater  than  10'^,  there  was  a  significant  amount  of  lag  and 
retardation  in  the  meniscus  position,  when  compared  to  an  inviscid  solution. 

The  effect  of  heating  the  capillary  tube  on  the  dynamic  response  was  also  smdied. 
Experimental  results  and  analytical  calculations  were  compared  for  a  Bond  munber  of 
1.47  and  a  capillary  number  of  3 .00x10'®.  Experimental  data  of  the  meniscus 

recession  and  advance  agreed  well  with  a  one-dimensional  equation  of  motion  accounting 
for  the  "dynamic"  contact  angle. 
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PREFACE 


Recent  technological  advances  in  such  areas  as  high  power  electronics  and  electric 
motor  design  have  provided  the  impetus  for  the  Air  Force  to  develop  more  electric 
technologies  in  high  performance  aircraft.  In  doing  so,  the  thermal  management  of 
electncal  components,  such  as  might  be  developed  for  actuation  systems,  may  play  a  key 
role  in  the  viability  of  a  more  electric  concept.  The  use  of  passive  thermal  management 
concepts  are  currendy  being  addressed  to  take  advantage  of  the  inherent  reliability  and 
simplicity  in  design.  Such  devices  may  incorporate  capillary  stractures  to  aid  heat 
transfer  and  have  been  generally  thought  of  as  being  urureliable  in  high  acceleration 
environments.  However,  litde  is  known  of  the  actual  performance  of  capillary  driven 
devices  in  a  transient  acceleration  environment  typical  of  high  performance  aircraft.  It 
is  the  objective  of  this  research  to  address  and  identify  fundamental  physical  concepts 
which  win  ultimately  aid  in  predicting  the  performance  and  design  parameters  of  capillary 
driven  heat  transfer  devices  subject  to  transient  accelerations. 

This  manuscript  is  the  compilation  of  three  papers  (chapters  I-DI)  for  publication 
in  the  Journal  of  Thermophysics  and  Heat  Transfer  and  the  9*  International  Heat  Pipe 
Conference.  As  a  result,  there  is  some  repetition  in  the  text  and  figures  to  maintain 
clarity  in  each  chapter.  The  reader  is  encouraged,  along  with  the  introduction  and  chapter 
IV  (summary  and  conclusions),  to  treat  chapters  I,  H,  or  HI  as  independent  investigations. 
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INTRODUCTION 


The  static  and  dynamic  behavior  of  two-  and  three-phase  interfaces  constitute  a 
range  of  problems  that  have  been  frequently  addressed  in  textbooks’"^  and  by 
investigators^"®  over  the  last  several  years.  These  interfaces,  which  may  be  referred  to  as 
contact  lines,  often  involve  the  displacement  of  one  immiscible  fluid  in  another  along  a 
solid  surface,  the  spreading  of  a  liquid  drop  on  a  solid,  the  displacement  of  gas  or  vapor 
by  a  liquid  in  a  ciqpdlary  tube,  etc.  There  are  many  processes  that  require  information 
with  regard  to  the  static  and  dynamic  behavior  of  contact  lines,  such  as  the  spreading  of 
viscous  liquids  like  paint,  adhesive  technologies,  lubrication  technologies,  oil  recovery, 
and  processes  utilizing  contact  lines  and  cj5)illarity  to  enhance  heat  and  mass  transfer. 
The  latter  relates  to  condensation  and  evaporation  processes  associated  with  a  wetting 
liquid  on  planar  surfaces  and  within  wicking  structures  consisting  of  porous  mpHia, 
grooved  stractures,  and  screens.  In  these  cases,  it  is  typical  to  have  length-scales 
associated  with  the  size  of  droplets,  grooves,  or  pores  on  the  order  of  several  millimeters 
or  less,  resulting  in  a  dominant  surface  tension  force  which  controls  the  shape  of  the 
fluid-fluid  interface  close  to  the  contact  line. 

Additional  difficulties  arise  when  devices  incorporating  these  processes  are 
subjected  to  elevated  steady  state  and  transient  acceleration  fields.  Elevated  steady  state 
acceleration  fields  resulting  in  forces  on  the  same  order  of  magnitude  or  greater  than 
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surface  tension  forces  (more  properly,  the  adhesion  forces  between  the  liquid  and  solid 
in  the  vicinity  of  the  near  contact  line  region)  may  affect  the  sh^  of  the  fluid  interface. 
Transient  acceleration  induced  forces  (both  transverse  and  axial),  in  addition  to  altering 
the  fluid  interface  shape,  will  cause  motion  of  the  contact  line  as  the  gravitational 
potential  associated  with  the  liquid  filled  pore  or  groove  changes.  In  such  instances,  the 
contact  line  velocity  is  functionally  dependent  on  these  transient  acceleration  induced 
forces  as  weU  as  surface  tension.  The  inclusion  of  these  elevated  steady  state  and 
transient  acceleration  fields  defines  a  class  of  problems  that  are  unique  in  the  study  of  the 
static  and  dynamic  behavior  of  contact  lines. 

If  the  problem  of  a  static  meniscus  in  a  capillary  tube  is  considered  (discounting 
the  near  contact  line  region  where  the  solid-liquid  intermolecular  forces  may  dominate  the 
force  field  at  the  fluid  interface),  the  equilibrium  equation 

V'd  +  pS^O  (1) 


must  be  satisfied.  A  traction  stress  vector  boimdary  condition  in  terms  of  surface  tension 
is  imposed  on  the  free-surface  interface.  If  this  interface  is  assumed  to  have  a  ^herical 
sht^ 


"(jf) 


n  = 


o 


(2) 


is  also  satisfied.  With  this  assumption,  the  "apparent"  contact  angle  at  the  contact  line 
is  also  required  as  a  boimdary  condition  to  determine  the  radius  of  curvature  and  therefore 
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the  interface  sh^e.  The  specification  of  an  apparent  contact  angle  is  intended  to 
macroscopically  describe  the  contribution  of  the  physics  of  the  extended  meniscus  region. 

However,  the  contact  angle  as  described  by  Young®  is  dependent  on  the  free 
energy  state  at  the  contact  line.  Yoimg  was  the  first  to  derive  an  expression  for  contact 
angle  using  a  force  balance  at  the  contact  line  and  represented  by  the  expression 

Oi^cosQ^  =  (3) 

This  shows  that  the  contact  angle  is  functionally  dependent  on  the  free  energy  state  at  the 
contact  line  and  may  be  used  as  an  indirect  determination  of  this  minimum  free  energy 
state.  Furthermore  the  difference  between  the  free  energy  of  the  liquid-solid  and  liquid- 
gas  to  that  of  the  solid-gas  is  defined  as  a  spreading  parameter 

‘^=  +  (4) 

or 

S  =  o^^(cose^-l)  (5) 

When  5 < 0,  the  solid-gas  interface  has  a  lower  free  energy  than  the  liquid-solid  and 
liquid-gas  and  therefore  the  free  energy  is  already  minimized  and  the  liquid  will  not 
spread. 

Although  static  contact  angles  have  been  reported,  early  experimenters  such  as 
Padday  used  direct  methods  of  photography  and  indirect  methods  using  measurements 
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from  the  vertical  position  of  the  contact  line  with  respect  to  the  xindisturbed  free  surfaces. 
Others  attempted  to  deduce  its  value  from  dynamic  experimental  measurements  using 
tangent  measurements  along  the  free-surface  at  a  point  removed  from  the  contact  line. 
These  investigations  treated  the  contact  line  as  a  region  of  zero  thickness  intersecting  the 
solid  wall  along  a  line  with  constant  surface  tension  and,  as  such,  were  unable  to  consider 
the  physics  governing  the  contact  line  region.  Hansen  and  Toong“  concluded  that 
classical  laws  of  fluid  motion  were  not  applicable  at  the  contact  line  and  that  observed 
contact  angles  may  not  be  the  appropriate  contact  angle  of  interest. 

Concus*^  investigated  the  effect  of  Bond  number  on  the  shape  of  this  interface  in 
a  right  circular  cylinder.  Several  contact  angles  were  specified  and  the  equilibrium 
equation  was  solved  for  the  interface  shape  as  a  function  of  contact  angle.  His  results 
showed  that  the  interface  shape  was  not  necessarily  spherical  in  nature,  but  was  strongly 
dependent  on  the  specified  contact  angle  and  corresponding  Bond  number. 

Several  different  approaches,  either  macroscopic  or  microscopic,  have  been  used 
to  describe  the  dynamics  associated  with  an  advancing  or  receding  uiiheated  meniscus. 
In  the  macroscopic  approach,  the  details  of  the  near  contact  line  region  are  ignored. 
Elliott  and  Riddiford*^  and  Philhps  and  Riddiford*^  investigated  the  advancing  and 
receding  contact  angle  for  a  growing  bubble  between  two  flat  plates  as  a  function  of 
interfacial  velocity.  They  found  a  definitive  lelationsh^  between  the  advancing  or 
receding  contact  angle  on  interfacial  velocity.  For  single  liquid  systems  spreading  on  a 
dry  surface,  the  advancing  contact  angle  was  found  to  be  independent  of  the  velocity  for 
interfacial  velocities  up  to  1  mm/min.  For  interfacial  velocities  greater  than  1  mm/min.. 
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the  advancing  contact  angle  increased  linearly  to  a  limiting  value.  Qualitative  molecular 
considerations  were  used  to  explain  the  velocity  dependence  of  advancing  contact  angles. 
Relaxation  of  the  contact  angle  from  a  dynamic  value  to  a  static  value,  when  the  fluid 
drive  system  was  removed,  was  described  as  a  molecular  reorientation  of  the  liquid 
molecules  at  the  contact  line. 

Hoffman*^  investigated  the  shape  of  an  advancing  interface  in  a  capillary  tube  for 
which  viscous  and  interfacial  forces  are  the  dominant  salient  features  defining  the  system. 
He  found  a  relationship  between  the  advancing  contact  angle  and  the  capillary  number 
(Ca  =  pV/o)  plus  a  shift  factor  determined  by  the  static  contact  angle. 

Recently,  Calvo  et  al.‘®  examined  the  advancing  meniscus  in  a  horizontal  capillary 
tube  attached  to  a  liquid  reservoir.  It  was  found  that  the  advance  of  the  meniscus  was 
controlled  by:  (i)  the  gravity  head  of  the  liquid  above  the  c^illary  tube;  and  (ii)  the  line 
force  associated  with  the  parent  contact  line  equal  to  (2o/r)  cos0  where  0  is  the 
apparent  contact  line  angle  for  dynamic  conditions.  Their  results  indicated  that  for 
capillary  numbers  less  that  10■^  0^  -  0^,  the  static  contact  angle.  With  increasing 
capillary  numbers  the  dynamic  contact  angle  increases  for  advancing  menisci  and 
decreases  for  receding  menisci.  The  macroscopic  predictions  of  the  advancing  rate  of  the 
meniscus  showed  good  agreement  with  experimental  H^ta, 

The  microscopic  approach  has  focused  on  the  details  of  the  thin  film  region  near 
the  contact  line.  Instead  of  prescribing  a  line  force,  this  ^proach  has  sought  to  determine 
the  total  dissipation  in  the  thin  film  where  the  velocity  gradients  will  be  maximum.  Thus, 
the  viscous  losses  in  the  thin  film  will  control  the  motion  of  the  meniscus.  Relative  to 
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these  microscopic  approaches,  de  Gennes  et  al.*^  Ngan  and  Dussan,  V.‘®,  and  Dussan,  V. 
et  al.  have  performed  many  of  the  most  recent  investigations  concerning  contact  line 
dynamic  behavior.  In  particular,  de  Gennes  et  al.®  considered  the  dynamic  behavior  of 
the  contact  line  associated  with  the  spreading  of  a  partially  wetting  fluid  on  an  inclined 
moving  plate.  They  distinguished  four  sub-regions,  as  shown  in  Fig.  1,  consisting  of: 
(i)  a  molecular  region  where  the  continuum  description  is  no  longer  valid;  (ii)  a  proximal 
region  dominated  by  the  long-range  van  der  Waals  forces;  (iii)  a  central  region  dominated 
by  viscous  and  capillary  forces;  and  (iv)  a  distal  region  dominated  by  viscous  and 
gravitational  body  forces.  Their  solution  approach  was  to  use  matched  as)miptotic 
expansions  to  describe  each  distinct  sub-region  such  that  the  boundary  conditions  of 
neighboring  sub-regions  were  satisfied.  The  problem  for  each  sub-region  was  formulated 
such  that  the  physics  was  correctly  defined  (i.e.,  proper  length-scales  and  ^ropriately 
dominant  salient  features).  They  solved  for  the  static  fluid  thickness  in  the  proximal  and 
central  region,  then  considered  the  contact  line  velocity  as  a  perturbation  and  determined 
a  first  order  correction  to  the  static  profile.  Finally  they  matched  this  profile  to  that  of 
the  distal  region  profile. 

Ngan  and  Dussan,  V.‘®  presented  an  investigation  whose  objective  was  to 
determine  the  appropriate  boundary-value  problem  in  an  outer  region  (analogous  to  de 
Gennes’  central  and  distal  regions)  excluding  the  inner  region  (analogous  to  de  Gennes’s 
proximal  and  central  regions)  in  the  immediate  vicrmty  of  the  moving  contact  line 
between  two  parallel  glass  plates.  The  formulation  of  this  boundary-value  problem  was 
based  on  the  predication  that  the  velocity  field  has  an  asymptotic  structure  near  the 
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moving  contact  line.  It  was  proposed  that  this  structure  contains  only  one  parameter 
made  up  of  terms  consisting  of  the  experimentally  measurable  quantities  of  apex  height 
and  intermediate  contact  angle,  0  at  r  =  J?  (slightly  removed  from  the  contact  line).  The 
boundary-value  problem  was  posed  such  that  its  validity  in  the  outer  region  was 
maintained  whale  a  slip  boundary  condition  was  applied  in  the  inner  region  along  the  solid 
surface.  Further  conditions  required  that  the  Bond  number,  capillary  number,  and 
Reynolds  number  approach  zero.  Their  investigation  dealt  primarily  with  the  extent  to 
which  this  formulation  accurately  described  the  shr^  of  the  fluid  interface  asyn^totically 
near  the  moving  contact  line  in  the  outer  region  as  r  -  0 ,  and  for  small  capillary  numbers, 
Ca  -*  0 .  The  validity  of  this  approach  was  verified  in  conjunction  with  experimental  data. 

Ngan  and  Dussan,  V.**  also  found  a  systematic  deviation  between  the  predicted 
and  experimental  observations  due  in  part  from  the  fact  that  the  theory  accounted  for  only 
the  two  lowest  orders  in  the  as3miptotic  e^ansion  in  the  capillary  number  and  hence 
indicated  that  the  capillary  number  used,  on  the  order  of  10'^,  may  have  been  too  large. 
A  subsequent  experimental  investigation  was  performed  by  Dussan,  V.  et  al.*’  to  verify 
these  hypotheses.  Experiments  were  performed  for  capillary  numbers  ranging  from 
2.8  X 10"^  to  8.3  X 10'^.  The  capillary  numbers  used  for  this  investigation  were  lower  by 
an  order  of  magnitude  than  previous  investigations  in  an  attempt  to  verify  that  the 
systematic  deviation  between  experimental  results  and  analytical  calculations  was  due  to 
the  neglect  of  higher-order  terms  in  the  capillary  number  expansions.  They  found  good 
agreement  between  the  experimental  results  and  the  theoretical  predictions  with  the 
systematic  deviations  appearing  only  at  the  higher  capillary  numbers. 
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With  the  addition  of  heat  to  the  meniscus  region,  the  dynamics  associated  with  an 
advancing  or  receding  meniscus  may  be  significantly  altered  from  that  of  the  unheated 
case.  In  addition  to  evaporation,  surface  tension  will  decrease  with  increased  temperature 
and  temperature  gradients  will  result  in  a  spatial  surface  tension  gradient  along  the 
meniscus  surface.  These  surface  tension  gradients  cause  interfacial  tangential  stresses  or 
thermocapillary  stresses  to  develop  at  the  meniscus  surface.  The  induced  stresses  due  to 
thermocapdlary  effects  have  the  effect  of  altering  the  contact  angle.  Modifying  the 
equation  by  Youngf,  eq.  (3),  to  account  for  the  thermocapillary  stress,  results  in 

cosQ^  =  (6) 


If  non-isothermal  conditions  exist,  thermoc^illary  stresses  at  the  interface  may  have  a 
net  effect  of  either  increasing  or  decreasing  the  contact  angle  depending  on  the  direction 
of  the  spatial  surface  tension  gradient. 

There  have  been  many  noteworthy  investigations  widi  regard  to  the  evaporating 
meniscus  with  and  without  the  inclusion  of  thermocqjfllary  effects.^®-"  Deijaguin  et  al.“ 
was  first  to  develop  a  theory  of  evaporation  from  a  ctq)illary  with  verification  from 
experimental  results.  Holm  and  Goplen^‘  developed  an  approximate  method  for  predicting 
the  surface  heat  transfer  coefficients  in  the  vicinity  of  the  extended  meniscus  (analogous 
to  de  Gennes’  proximal  and  central  regions)  in  a  capillary  groove  stracture.  They 
concluded  that  heat  transfer  through  the  extended  meniscus  accoimted  for  approximately 
80%  of  the  total  heat  dissipation  from  the  wall.  Moosman  and  Homsy^  used  perturbation 
theory  to  model  a  horizontal  non-isothermal  extended  meniscus.  They  showed  that  there 
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is  a  potential  for  large  heat  fluxes  to  occur  near  the  contact  line  with  a  subsequent 
thinning  of  this  region  and  an  increase  in  the  apparent  contact  angle.  However,  none  of 
these  investigators  considered  thermocapiUary  effects. 

Investigations  taking  into  account  thermocapillary  effects  have  been  primarily 
mathematical  in  nature  without  experimental  verification.  Sen  and  Davis^  analytically 
demonstrated  thermocapiUary  flows  in  a  two-dimensional  slot  induced  as  a  result  of  an 
applied  temperature  gradient  to  the  free  surface.  They  assumed  either  a  fixed  contact  line 
or  a  fixed  contact  angle  as  a  basis  for  their  analytical  solutions.  Swanson  and  Herdt^‘‘ 
developed  a  mathematical  model  of  an  evaporating  meniscus  in  a  capUlary  tube  talcing 
into  account  thermocapiUary  effects  for  smaU  heat  inputs.  They  observed  no  change  in 
the  apparent  contact  angle  and  noted  that  the  meniscus  profile  remained  spherical. 
HaUinan  et  al.^  considered  the  same  geometry  as  Swanson  and  Herdt  also  allowing  for 
non-isothermal  extended  meniscus.  Results  were  comparable  with  other  investigators  and 
they  identify  conditions  necessary  for  thermocapiUary  stresses  to  have  an  effect  on  the 
flow  into  the  extended  meniscus  region. 

Pratt  and  Hallinan“  present  an  experimental  and  analytical  approach  on  the 
thermocapiUary  effects  on  a  heated  meniscus  in  a  c^iUaiy  tube.  They  experimentaUy 
demonstrate  their  claim  that  thermocapiUaiy  stresses  induce  a  meniscus  recession  and 
present  the  critical  A  T  or  Marangoni  number  required  for  the  meniscus  to  recede.  The 
critical  A  T  was  analyticaUy  shown  to  have  the  effect  of  increasing  the  apparent  contact 
angle  resulting  in  a  meniscus  recession. 

To  date,  these  formulations  have  been  limited  to  situations  where  the  Bond  number 
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is  small  (e.g.,  inside  of  capillary  tubes)  or  where  it  is  infinite  (such  as  exists  for  draining 
films^’).  While  acceleration  induced  forces  are  considered  conservative  and  long-range 
in  nature,  they  affect  the  free  surface  shape  and  the  slope  for  which  solutions  to  the  sub- 
regions  as  described  by  de  Geimes  et  al.‘’  must  be  asymptotically  matched.  Furthermore, 
it  is  not  known  what  potential  effects  these  forces  may  have  on  the  precursor  film 
thickness  and  stability  in  both  static  and  dynamic  conditions. 

In  general,  prior  investigations  with  regard  to  the  dynamic  behavior  of  the  contact 
line  have  specified  a  contact  line  velocity  consistent  with  a  moving  solid  boundary. 
Transient  acceleration  induced  forces  wiU  induce  a  contact  line  velocity  posing  a  problem 
such  that  the  contact  line  d3aiamics  is  dependent  only  on  the  temporal  acceleration  fields 
and  not  on  a  moving  solid  boimdary.  There  has  been  limited  research  with  regard  to  the 
effects  of  acceleration  induced  forces.  Most  of  these  have  addressed  the  equilibrium  fiee- 
suiface  shape. This  investigation  addresses  the  dynamic  behavior  of  a  unheated  and 
heated  meniscus  and  contact  line  region  in  a  capiSidiry  tube  subjected  to  a  transient 
acceleration  field. 
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CHAPTER  I 


DYNAMIC  BEHAVIOR  OF  AN  UNHEATED  MENISCUS 
SUBJECTED  TO  A  TRANSIENT  ACCELERATION  HELD 

Abstract 

This  study  experimentally  and  analytically  investigates  the  effects  of  acceleration 
transients  on  the  dynamic  response  of  an  unheated  meniscus  in  a  capUlary  tube.  A 
capillary  tube  partially  filled  with  either  water  or  ethyl  alcohol  was  mounted  on  a 
centrifuge  to  observe  the  dynamic  response  of  the  meniscus  subject  to  an  acceleration 
transient.  It  was  found  that  for  a  transverse  acceleration  component  of  less  than  2g, 
experimental  data  of  the  meniscus  advance  or  recession  agreed  well  with  a  one¬ 
dimensional  formulation  of  the  equation  of  motion  which  assumed  a  spherical  shtqred 
meniscus.  However,  for  transverse  accelerations  greater  than  2g,  there  was  an  increase 
in  the  observed  meniscus  velocity  both  upon  recession  and  advancement  resulting  in  a 
lower  meniscus  position  when  compared  to  the  analytically  predicted  position.  Water 
exhibited  a  signiricant  amount  of  "sticking"  and  hysteresis  which  appeared  to  intensify 
with  increasing  transverse  acceleration. 
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Experimental 


Acceleration  Field 

Transient  acceleration  fields  were  generated  using  a  2.4-m  diameter  centrifuge 
table  (Fig.  1.1)  rotating  with  a  time  variant  angular  velocity.  Angular  velocity  transients 
were  generated  at  a  specified  cyclic  frequency  of  0.0015  Hz(t£  =  I/O. 0015  s)  using 
a  signal  or  waveform  generator  as  a  control  input  to  the  centrifuge.  The  angular  velocity 
was  increased  linearly  to  a  peak  value  at  the  midpoint  of  the  cycle  and  subsequentiy 
decreased  as  shown  in  Fig.  1.2a.  A  tri-axis  accelerometer  (Columbia  Research 
Laboratories,  Inc.)  was  used  to  monitor  the  time  variant  acceleration  components  in  a 
cartesian  reference  frame  affixed  to  the  centrifuge  table.  Figure  1.2b  shows  a  typical 
accelerometer  output  to  the  experimentally  generated  acceleration  field  for  each  of  the 
three  coordinate  directions.  For  this  investigation,  the  peak  angular  velocity  was  specified 
such  that  a  4.2-g  ±  0.1-g  peak  radial  acceleration  was  generated  (gr^jj,  =  41.20  m/s^). 

Capillary  Test  Cell 

A  sealed  test  cell  containing  a  glass  capillary  tube  and  reservoir  was  positioned 
on  the  centrifuge  as  shown  in  Fig.  1.3.  The  test  cell  was  mounted  to  a  motorized  optics 
rotation  stage  (Newport)  such  that  the  capillary  tube  was  allowed  to  pivot  about  its  center 
of  mass.  This  pivot  point  was  displaced  vertically  from  the  centrifuge  surface  at  a  fixed 

radial  location,  r^,  as  shown  in  Fig.  1.4a.  Angles  of  inclination  of  the  capillary  tube 

were  determined  from  the  output  of  a  calibrated  angular  displacement  transducer  (Trans- 
Tek,  Inc.)  accurate  to  within  ±0.2  degrees. 
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Experiments  were  conducted  by  fixing  the  capillary  tube  at  a  specified  inclination 
angle  while  subjecting  it  to  the  aforementioned  transient  acceleration  field.  Increasing  the 
inclination  angle  elevates  the  transverse  acceleration  component  referenced  to  the  capillary 
tube.  This  allowed  the  dynamic  response  of  the  meniscus  subjected  to  an  mcreasing 
transverse  acceleration  component  to  be  observed  experimentally.  The  dynamic  response 
of  the  meniscus  to  the  transient  acceleration  field  was  observed  using  an  8-mm  CCD 
camera  (Sony)  mounted  adjacent  to  the  test  cell  as  shown  in  Fig.  1.3b.  The  meniscus 
height  relative  to  the  reservoir  meniscus,  t]  ,  was  measured  as  a  function  of  time  to  within 
±1.0  mm.  An  equivalent  contact  angle,  0,  was  calculated  from  the  experimentally 
determined  static  wicking  height,  hg,  for  conditions  when  the  transverse  acceleration 
component  was  zero  ( h^p  g  -  (2o/a)cos0). 

Water  (reverse  osmosis)  and  ethyl  alcohol  were  used  as  the  test  fluids  to 
demonstrate  the  importance  of  the  equilibrium  contact  angle  on  the  advance  and  recession 
of  the  meniscus  due  to  a  transient  acceleration  field.  Properties  of  these  fluids  were 
assumed  to  be  that  of  the  bulk  reservoir  fluid  corresponding  to  a  experimental  temperature 
of  3(rc. 

The  test  cell  and  capillary  tube  were  cleaned  using  potassium  hydroxide  and  rinsed 
with  distilled  water.  The  test  cell  was  filled  and  subsequently  sealed  resulting  in  a  fill 
that  had  a  combination  of  air,  fluid  v^or,  and  liquid. 
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Analytical  Fonnulation 


The  goals  of  the  analytical  formulation  were  to;  (i)  mathematically  describe  the 
transient  nature  of  the  acceleration  field  referenced  to  the  capillary  tube  in  the  meniscus 
region;  (ii)  mathematically  describe  the  dynamic  behavior  of  the  meniscus  as  a  result  of 
the  temporal  acceleration  induced  forces  using  a  simplified  one-dimensional  model  with 
a  constant  line  force;  and  (iii)  determine  the  relative  importance  of  the  transverse 
acceleration  component,  Bond  number,  and  capillary  number  on  the  meniscus  motion  by 
comparing  experimental  results  and  analytical  calculations. 

Acceleration  Vector 

Acceleration  measurements  were  obtained  with  regard  to  a  three-dimensional 
cartesian  non-inertial  reference  frame  affixed  to  the  centrifuge.  The  acceleration  field  was 
then  referenced  to  fixed  locations  on  the  capillary  tube  and  subsequently  deconq}osed  into 
transverse  and  axial  acceleration  components  relative  to  the  capillary  tube  as  shown  in 
Fig.  1.4b. 

The  transient  acceleration  vector  in  a  cartesian  coordinate  reference  frame  from 
the  viewpoint  of  an  observer  on  the  centrifuge  is  of  the  form 

=  b^§^+b^§j^  (1) 

where  for  a  fixed  radial  location,  r  the  following  acceleration  components  are  realized: 

radial  component,  b^^(t)  =w^(t)rj., 
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tangential  component, jbj.(  t)  = 
vertical  component,  hy{t)  =  -g. 


After  an  appropriate  transformation  to  the  capillary  tube  reference  frame,  the  acceleration 
vector  can  be  decomposed  into  an  axial  component,  and  a  transverse  component,  2)^, 
which  results  in 


b  =  [-Ci)^(ro+  (z^-z)  costf))  cos<l>+  ( -gr)  sin<j)] + 
[-0)^  (ro+ (Zi-z)  cos4>)  sirwj)- (-gr)  cos<j)]^  +  (ro+ (Zj^- z)  cos<J)) 


i/j 


(2) 


Here,  the  transverse  component,  jb^j,  is  a  magnitude  with  no  reference  to  direction  due  to 
the  axisymmetric  nature  of  the  capillary  tube. 

Equation  of  Motion 

A  simple  analytical  model  is  described  which  predicts  the  motion  of  the  meniscus 
subjected  to  transient  accelerations.  The  analytical  formulation  of  the  equation  of  motion 
was  simplified  assuming  bulk  flow  only  in  the  axial  direction  of  the  capillary  tube. 
Effectively  the  liquid  column  in  the  capillary  tube  is  assumed  to  undergo  a  slug  flow. 
From  the  conservation  of  momentum, 

Ef  =  (3) 
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the  surface  and  body  forces  in  the  axial  direction  combine  resulting  in  the  axial  force 

(4) 

The  surface  force,  ,  is  formulated  by  accounting  for  the  cumulative  effect  of 

solid-liquid  intermolecular  forces  in  the  near  contact  line  region  and  the  wall  shear  stress 
associated  with  Poiseuille  flow  in  the  liquid  column.  The  contact  line  force  is  defined 
assuming  a  spherical  meniscus  and  a  constant  equivalent  or  apparent  contact  angle  such 
that 

Fg^  =  o(27ta)cos0  +  T^(27ta)ri.  (5) 


Note  that  such  an  assumption  implies  a  non-changing  apparent  contact  angle  under 
dynamic  conditions.  A  great  deal  of  work  has  established  a  relationship  between  the 
dynamic  contact  angle,  0  and  the  capillary  number  (for  low  Bond  numbers)  as 
summarized  by  Hoffinan*^.  These  results  have  shown  little  change  in  the  contact  angle 
for  capillary  numbers  less  than  10'^. 

The  shear  stress  at  the  wall  is  defined  using  the  Darcy  friction  factor  assuming 
laminar  tube  flow  where 

^  64  (6) 

Re[  8  ) 

This  results  in  the  following  surface  force  formulation  while  using  a  slug  flow 
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approximation  to  describe  the  flow  of  liquid  in  the  column,  i.e.,  =  dr\/dt. 


F  =  0(2X3)0030-87:^11 


dt 


(7) 


The  total  body  force, 

(2x)  rdrdz. 


(8) 


is  formulated  by  integrating  over  the  volume  of  the  liquid  coltimn  the  axial  acceleration 
component,  from  eq.  (2)  as  defined  by 

b^  =  -(o(t)^  cost})  (ro+  cos4>)-srsin<|).  (8) 


The  transient  angular  velocity  is  of  the  form 


(0  ( t)  =  Bt, 

=  B(  tf-  t)  , 


tf 

2 


tf 


^t^tf 


(10) 


where 


B  = 


^  ^pk 


itf) 


2  r' 
O 


1 

2 


(11) 


Assuming  that  the  axial  fluid  velocity  is  also  defined  by  =  dq/dt  (due  to  the  slug 
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flow  assumption)  and  the  mass  in  the  capillary  tube  asm=  p(7ia^)Ti,the  relation  for 
the  time  rate  of  change  in  momentum  is  defined  as 


dimvj 

dt~ 


P 


(Tza^) 


(12) 


Combining  eqs.(7),  (8),  and  (12),  the  momentum  equation  becomes 


■  2OCOS0 

pa 

pa2  '  dt\  [Jo  ^  J 

['  dt2  \  dtj  \ 

(13) 


Integrating  the  axial  acceleration  term,  h^,  and  rearranging,  the  momentum  equation 
becomes 


2ocos6 

pa 


r  dt^  \dcj\ 

[pa^  dtJ 

cos<|)|(Zi- -^)  cos4)  +  ro|  +  gsin<|)  , 


(14) 


A  non-dimensional  form  of  the  momentum  equation  can  now  be  derived  using  the 
following  dimensionless  parameters 


C  =  0  ^  C  ^  1 

ho 

©  =  0:^  ©  ^  1, 


(15) 


with  the  following  initial  conditions 
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C  =  1,  -^  =  0,0  =  0. 


A  new  angular  velocity,  Q ,  is  formed  by  factoring  out  of  the  coefficient  B  (defined 
by  eq.  (11))  in  eq.  (10)  and  casting  it  in  terms  of  the  non-dimensional  time,  0,  where 


4  \i 


£2(0)  =  ^ 

[roj 


0, 


(1-0), 


0  :^0:S  -5^ 
2 

— 

2 


(16) 


The  non-dimensional  form  of  the  momentum  equation  in  terms  of  a  transient 
dimensionless  meniscus  position  becomes 


-^sin<b 


a 

■  ReCa' 

2Bo 

d02  \  d0  /  J  [  So 

L 

£2  (0)  ^  cost}) 


(17) 


The  first  and  second  terms  of  eq.  (17)  represent  the  inertial  and  viscous  effects, 
respectively,  with  the  coefficients  formed  as  the  products  of  the  dimensionless  parameters 
of  Reynolds,  capillary,  and  Bond  numbers.  The  Reynolds  and  capillary  numbers  are 
referenced  to  a  characteristic  velocity,  v^,  which  is  the  maximum  attainable  velocity  over 
one-half  of  the  acceleration  cycle  period  assuming  q  =  0  at  The  Bond  number  is 
referenced  to  the  peak  radial  acceleration  component, 

For  small  inertial  and  viscous  effects,  the  momentum  equation  can  be  reduced  to 
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the  following  quadratic  form  for  the  transient  dimensionless  meniscus  position. 


Q  (0)  ^  cos<l>  (z,cos<j)  +  r'o)  +-S-  sin<j) 
_ 9pk 

£2  (©)  ^  —  cos^4) 

2 


c  + 


— ^sind) 

_ _  _  I 

Q  (0)  ^  —  cos^<|) 

2  ; 

(18) 


Results 

Equations  (17)  and  (18)  were  solved  for  the  dimensionless  meniscus  position  using 
a  central-difference  scheme  and  as  a  quadratic  respectively.  These  solutions  were 
obtained  as  a  result  of  specifying  a  static  wicking  height,  hg,  based  on  the  experimentally 
determined  equivalent  contact  angle,  0 .  Solutions  to  eqs.  (17)  and  (18)  were  identical  for 
the  experimental  parameter  range  specified  in  this  investigation.  This  supports  the  use 
of  the  closed  form  solution  of  eq.  (18)  in  instances  where  inertial  and  viscous  effects  can 
be  assumed  to  be  negligible. 

Experiments  were  performed  and  the  results  compared  with  solutions  to  the 
analytical  model.  Uncertainty  in  the  experimentally  obtained  dimensionless  meniscus 
position,  C»  and  time,  0,  was  on  the  order  of  C,  ©  ±  0.02  based  on  the  accuracy  and 
resolution  of  the  experimentally  measured  quantities.  In  general,  the  comparison  between 
the  model  and  the  experiment  is  surprisingly  good,  given  the  simplicity  of  the  model. 
However,  a  number  of  interesting  deviations  from  the  model  help  to  provide  insight  into 
the  physical  phenomena  associated  with  the  advancing/receding  meniscus. 

Figure  1.5  compares  the  experimentally  obtained  dimensionless  meniscus  position 
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to  the  analytical  results  for  water  at  varying  angles  of  inclination.  Also  shown  are  the 
axial  and  transverse  acceleration  components.  The  equivalent  contact  angle  was  22.2® 
with  an  experimental  Bond  number  of  2.05.  Experimental  parameters  ranged  from 
5.9X10"''  ^  Ca  ^  1.5xl0-«,  0.11  ^  J?e  ^  0.28,  and 
2.3x10”^  ^  e-lh^  ^  5.7x10”^  as  the  inclination  angle  varied  from  20°  to  60°. 

There  was  a  significant  amount  of  "sticking"  of  the  advancing  meniscus  for  water 
as  evidenced  by  the  discontinuous  nature  of  the  meniscus  position  curves,  Fig.  1.5.  There 
also  was  a  permanent  hysteresis  in  the  meniscus  position  at  the  completion  of  the  test 
cycle  resulting  in  a  lower  wicking  height.  This  "sticking"  of  the  meniscus  appeared  to 
become  more  significant  with  an  increased  transverse  acceleration  component,  as 
demonstrated  in  Fig.  1.5e.  For  a  peak  transverse  acceleration  component  greater  than 
2.5g  the  experimental  results  showed  an  increase  in  the  meniscus  velocity  for  both  the 
receding  and  advancing  meniscus  from  that  predicted  by  the  analytical  results.  Figs.  1.5d 
and  5e.  At  the  midpoint  of  the  acceleration  cycle  the  meniscus  position  was  also  lower 
in  Figs.  1.5d  and  1.5e  than  that  predicted  analytically. 

Also,  for  all  cases  using  water  as  the  test  fluid  the  initial  meniscus  recession  speed 
tended  to  be  under-predicted.  Since  the  contact  line  force  in  the  predictions  is  controlling 
the  receding  and  advancing  velocities  it  is  ppaient  that  it  is  slightly  over-predicting  the 
actual  contact  line  dissipation.  Using  an  equivalent  contact  angle  to  describe  the  line 
force  at  the  contact  line  neglects  the  fact  that,  as  the  meniscus  recedes,  some  fluid  is  left 
behind  rather  that  being  dragged  over  the  tube  wall.^^  Thus,  the  line  force  is  physically 
less  than  predicted  by  the  model.  Conversely,  for  the  advancing  contact  line,  since  it  is 
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presumed  that  a  film  of  liquid  is  already  present  on  the  inside  diameter  of  the  tube,  it  is 
easier  for  the  meniscus  to  advance  relative  to  that  predicted  for  an  advancing  meniscus 
over  a  dry  surface  (such  as  is  associated  with  using  the  equivalent  contact  angle). 

Figure  1.6  compares  the  experimentally  obtained  dimensionless  meniscus  position 
to  the  analytical  results  for  ethyl  alcohol  at  varying  angles  of  inclination.  As  with  Figure 
1.5,  also  shown  are  the  axial  and  transverse  acceleration  components.  The  equivalent 
contact  angle  was  16.7°  with  an  ejqjerimental  Bond  number  of  1.47.  Experimental 
parameters  ranged  from  1.7x10'®  ^  Ca  ^  3.6x10'®, 

3.0x10'2  <;  i?e  ^  6  .lxlO'2,  and  1.9x10'^  ^  a/fi^  ^  4.0x10'^. 

The  experimentally  obtained  dimensionless  meniscus  position  for  ethyl  alcohol 
agrees  well  with  the  analytical  results  to  within  experimental  error  for  a  peak  transverse 
acceleration  component  less  than  2g.  No  "sticking”  of  the  meniscus  was  observed  and 
there  was  no  observable  hysteresis.  As  shown  in  Figs.  1.6  c-e,  there  was  an  observable 
increase  in  the  meniscus  velocity  relative  to  the  model  for  bodi  the  receding  and 
advancing  meniscus  for  a  peak  transverse  acceleration  component  greater  than  2g.  This 
increase  in  meniscus  velocity  resulted  in  a  meniscus  position  lower  than  that  predicted 
analytically. 

The  effects  of  capillary  number  for  a  fixed  Bond  number  and  inclination  angle  on 
the  dimensionless  meniscus  position  are  shown  in  Fig.  1.7.  A  comparison  of  the 
analytical  results,  using  eq.  (17),  for  a  Bond  number  of  0.42  and  capiUaiy  numbers  of 
lO"^,  10'^,  and  lO"^  at  a  20°  inclination  angle  demonstrates  that,  for  capillary  numbers  less 
than  10'^,  viscous  losses  outside  of  the  contact  line  had  only  a  small  effect  on  the 
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meniscus  position.  However,  for  capillary  numbers  on  the  order  of  lO'^  and  greater  there 
is  a  significant  lag  in  the  meniscus  position  due  to  viscous  effects.  C^illaiy  numbers  less 
than  10'^  and  Bond  numbers  ranging  fi-om  0.42  to  2.1  were  predicted  to  have  no  effect 
on  the  dynamic  response  of  the  meniscus.  Analytical  results,  using  eq.  (17),  for  a 
capillary  number  of  10"*  and  Bond  numbers  of  0.42,  0.84,  and  1.26  at  angles  of 
inclination  ranging  from  20”  to  60”  are  con^ared  in  Fig.  1.8.  Figure  1.8a  shows  the 
dimensionless  meniscus  position  at  varying  angles  of  inclination  for  a  fixed  Bond  number 
of  0.42  and  a  capillary  number  of  10"*.  For  decreasing  angles  of  inclination,  viscous 
effects  outside  the  contact  line  region  become  significant  resulting  in  the  characteristic  lag 
in  meniscus  position.  These  analytical  results  do  not  take  into  account  the  significance 
of  the  transverse  acceleration  component  at  the  higher  angles  of  inclination  and  to  the 
extent  at  which  this  transverse  component  affects  the  d5miamic  response  of  the  meniscus. 
Figure  1.8b  compares  the  dimensionless  meniscus  position  for  a  varying  Bond  number 
witii  that  of  the  inviscid  solution,  using  eq.  (18),  neglecting  inertial  effects  at  a  20” 
inclination  angle.  Here,  a  decrease  in  Bond  number  resulted  in  the  increased  importance 
of  viscous  effects  yielding  a  lag  in  the  meniscus  position  when  compared  to  the  low 
inertia  inviscid  solution.  For  Bond  numbers  less  than  1.3  and  capillary  numbers  greater 
than  10"*,  there  were  sigitificant  viscous  effects  resulting  in  a  lag  in  the  meniscus  position. 
However,  these  results  assume  a  constant  contact  line  force,  as  a  result  of  specifying  a 
constant  equivalent  contact  angle,  and  does  not  account  for  the  functional  dependence  of 
the  contact  angle  on  capillary  number  as  discussed  by  Hoffman*^.  However,  for 
experimental  and  analytical  results  with  capillary  numbers  less  than  10'^,  this  dependence 
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can  be  assumed  to  be  negligible. 

The  magnitude  of  the  transverse  acceleration  component  had  a  significant  effect 
on  the  dynamic  response  of  the  meniscus.  This  is  in  part  attributed  to  the  use  of  the 
constant  equivalent  contact  angle  to  define  a  constant  contact  line  force  in  the  analytical 
formulation.  While  the  meniscus  shape  may  remain  relatively  spherical  in  nature  with 
an  increased  axial  acceleration  component  for  the  relatively  low  Bond  numbers 
considered,  the  transverse  acceleration  component  wiU  serve  to  reorient  the  meniscus  in 
a  transverse  direction  and  distend  it  in  the  axial  direction.  This  variation  in  the  meniscus 
shape  and  position  in  the  capillary  mbe  will  result  in  an  increased  contact  line  length 
which  will  tend  to  be  elliptical  rather  than  spherical,  as  shown  in  Fig.  1.9.  Therefore,  one 
might  expect  an  increase  in  the  net  contact  line  force  due  to  the  increased  length  of  the 
contact  line.  Theoretically,  the  predicted  recession  and  advancing  rates  would  be  differ 
from  those  shown,  actually  moving  in  a  direction  opposite  of  that  seen  experimentally  in 
the  water  and  ethyl  alcohol  tests.  Such  observations  point  to  the  need  for  a  closer 
examination  of  the  near  contact  line  region  fluid  physics  with  regard  to  the  transverse 
acceleration  component  to  completely  describe  the  motion  of  the  meniscus.  Therefore, 
use  of  an  equivalent  contact  angle  to  define  the  contact  line  force  is  an  approximation  at 
best.  Instead,  the  contact  angle,  rather  than  being  constant,  is  a  function  of  the  transverse 
acceleration  component  and  angular  position  on  the  c^iUary  tube  wall. 
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a) 


Fig.  1.2  Plot  of  the  centrifuge  test  cycle  (0.0015  Hz):  a)  angular  velocity, 
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acceleration, 


b) 


Fig.  1.2(contd)  b)  resulting  acceleration  field. 
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Fig.  1.3  Experimental  apparatus:  a)  sealed  test  ceU 
and  b)  schematic  of  experimental  apparatus  as  mounted  on  the  centrifuge. 


a) 


Fig.  1.4  Capillaiy  tube  orientation:  a)  pivot  point  location,  directional  coordinates  and 
b)  coordinate  systems  as  referenced  to  the  centrifuge  accelerometer  and  capillary  tube. 
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Fig.  1.5  Experimental  data  compared  to  analytical  results  for  water  at  Bo  = 
a)  6  =  20°,  ca  =  1.50x10"®,  Re  =  0.278,  a/h„  =  2.25x10 


Fig.  L5(contd) 

b)  <|)  =  30®,  Ca  =  1.03x10'®,  Re  =  0.191,  a/hg  =  3.30x10 


acceleration, 


Fig.  1.5(contd) 

c)  A  =  40®,  Ca  =  7 .99x10-’,  Re  =  0.148,  a/l2„  =  4.24x10 


Fig.  1.5(contd) 

d)  <J)  =  50®,  Ca  =  6 .70xl0‘'^,  Re  =  0.124,  a/h.  =  5.05x10"^ 


acceleration, 
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Fig.  1.5(contd) 

e)<j)  =  60°,  Ca  =  5.93x10-’,  Re  =  O.llO,  a/h^  =  5.71x10 


acceleration, 


Fig.  1.6(contd) 

b)  <j)  =  30°,  ca  =  3.00x10'®,  Re  =  5.13x10-2,  a/h^  =  2.28x10 


I  f  ry 


acceleration, 


Fig.  1.6(contd) 

d)  <b  =  50°,  ca  =  1.96x10-®,  Re  =  3.35x10-2,  a//2„  =  3.50x10 


acceleration 


Fig.  1.7  Analytical  results  for  <j)  =  20°,  Bo  =  0.42  and  varying  ct5)iUary  number. 


CHAPTER  n 


DYNAMIC  BEHAVIOR  OF  AN  UNHEATED  MENISCUS 
SUBJECTED  TO  A  TRANSIENT  ACCELERATION  FIELD; 

EFFECT  OF  A  LOW  BOND  NUMBER 

Abstract 

This  study  experimentally  and  analytically  investigates  the  effects  of  acceleration 
transients  on  the  dynamic  response  of  an  unheated  meniscus  in  a  capillary  tube.  A 
capillary  tube  partially  filled  with  ethyl  alcohol  was  mounted  on  a  centrifuge  to  observe 
the  dynamic  response  of  the  meniscus  subject  to  single  and  multiple  cycle  acceleration 
transients.  Experimental  data  of  the  meniscus  recession  and  advance  agreed  well  with  a 
one-dimensional  equation  of  motion  accounting  for  changes  in  the  receding  and  advancing 
contact  angle  assuming  a  spherically  shaped  meniscus.  Systems  typified  by  a  low  Bond 
number,  as  distinguished  by  a  small  capillary  tube  diameter,  tended  to  be  insensitive  to 
transverse  acceleration  components.  It  was  found  that  the  combined  effects  of  a  Bond 
number.  Bo,  less  than  1.5  (distinguished  by  a  small  csqjiUary  tube  diameter),  and  capillary 
number,  Ca,  greater  than  lO'®  resulted  in  a  retardation  of  the  meniscus  position,  relative 
to  the  inviscid  solution,  due  to  acceleration  transients.  An  additional  retardation  in  the 
advancing  meniscus  position,  indicated  by  experimental  results,  appeared  to  be  due  to  an 
increase  in  the  advancing  contact  angle  with  increasing  capillary  number. 
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Experimental 


Acceleration  Field 

Transient  acceleration  fields  were  generated  using  a  2.4-m  diameter  centrifuge 
table  (Fig.  2.1)  rotating  with  a  time  variant  angular  velocity.  Angular  velocity  transients 
were  generated  at  a  specified  cyclic  frequency,  0. 0015  Hz  ^  f  .  10  Hz,  using  a 
signal  or  waveform  generator  as  a  control  input  to  the  centrifuge.  The  angular  velocity 
was  increased  linearly  to  a  peak  value  at  the  midpoint  of  the  cycle  and  subsequently 
decreased  as  shown  in  Fig.  2.2  for  single  and  multiple  cycles.  A  tri-axis  accelerometer 
(Coliunbia  Research  Laboratories,  Inc.)  was  used  to  monitor  the  time  variant  acceleration 
components  in  a  cartesian  reference  frame  affixed  to  the  centrifuge  table.  Figure  2.3 
shows  typical  accelerometer  outputs  to  the  experimentally  generated  acceleration  fields 
at  varying  fijequencies,  f,  for  single  and  multq)le  cycles  in  each  of  the  three  coordinate 
directions.  For  this  investigation,  the  peak  angular  velocity  was  specified  such  that  a  4.2- 
g  ±  0.1 -g  peak  radial  acceleration  was  generated  =  41.20  m/s^). 

C^iUaiy  Test  Cell 

A  sealed  test  cell  containing  a  glass  capillary  tube  and  reservoir  was  positioned 
on  the  centrifuge  as  shown  in  Fig.  2.4.  The  test  cell  was  mounted  to  a  motorized  optics 
rotation  stage  (Newport)  such  that  the  capillary  tube  was  allowed  to  pivot  about  its  center 
of  mass.  This  pivot  point  was  displaced  vertically  from  the  centrifuge  surface  at  a  fixed 

radial  location,  r^,  as  shown  in  Fig.  2.5a.  Angles  of  inclination  of  the  capillary  tube 

were  determined  from  the  output  of  a  calibrated  angular  displacement  transducer  (Trans- 
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Tek,  Inc.)  accurate  to  within  ±0.2  degrees. 

Experiments  were  conducted  by  fixing  the  capillary  tube  at  a  specified  inclination 
angle  while  subjecting  it  to  the  transient  acceleration  field.  Increasing  the  inclination 
angle  elevates  the  transverse  acceleration  component  referenced  to  the  capillary  tube. 
This  allowed  the  dynamic  response  of  the  meniscus  subjected  to  an  increasing  transverse 
acceleration  component  to  be  observed  experimentally.  The  dynamic  response  of  the 
meniscus  to  the  transient  acceleration  field  was  observed  using  an  8-mm  format  CCD 
camera  (Sony,  30  frames/s)  mounted  adjacent  to  the  test  cell  as  shown  in  Fig.  2.4b. 
Video  images  were  used  to  provide  a  record  of  the  meniscus  height  relative  to  the 
reservoir  meniscus,  ti  ,  as  a  function  of  time  to  within  ±1.0  mm.  An  equivalent  static 
contact  angle,  0^,  was  calculated  from  the  e:q)erimentaUy  determined  static  wicking 
height,  h^,  for  conditions  when  the  transverse  acceleration  component  was  zero 
{h^pg  -  (20/3)00365). 

Ethyl  alcohol  was  used  as  the  test  fluid  to  demonstrate  the  importance  of  the 
eqtiilibrium  contact  angle  on  the  advance  and  recession  of  the  meniscus  due  to  a  transient 
acceleration  field  for  low  Bond  numbers.  Two  capillary  tubes,  1.0  mm  and  0.5  mm  in 
diameter,  were  used  in  order  to  vary  the  experimental  Bond  number.  Properties  of  the 
ethyl  alcohol  were  assumed  to  be  that  of  the  bulk  fluid  corresponding  to  a  experimental 
temperature  of  30"  C. 

The  test  cell  and  capillary  tube  were  cleaned  using  potassiiun  hydroxide  and 
thoroughly  rinsed  with  ethyl  alcohol.  The  test  cell  was  filled  and  subsequently  sealed 
resulting  m  a  fill  that  had  a  combination  of  air,  fluid  vapor,  and  liquid. 
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Analytical  Formulation 


The  goals  of  the  analytical  fonnulation  were  to  expand  the  mathematical 
formulation  described  in  chapter  I  to:  (i)  mathematically  describe  the  dynamic  behavior 
of  the  meniscus  as  a  result  of  the  temporal  acceleration  induced  forces  using  a  simplified 
one-dimensional  model  with  a  constant  line  force  accounting  for  changes  in  the  recession 
and  advancing  contact  angles;  and  (ii)  determine  the  relative  importance  of  the  transverse 
acceleration  component  and  capillary  number  for  low  Bond  numbers  on  the  meniscus 
motion  by  comparing  experimental  results  and  analytical  calculations. 


Acceleration  Vector 

Acceleration  measurements  were  obtained  with  regard  to  a  three-dimensional 
cartesian  non-inertial  reference  frame  affixed  to  the  centrifuge  as  described  in  chapter  I. 

The  acceleration  field  was  then  referenced  to  fixed  locations  on  the  capUlaiy  tube  and 
subsequently  decomposed  into  transverse  and  axial  acceleration  components  relative  to 
the  capillary  tube  as  shown  in  Fig.  2.5b. 

The  resulting  transient  acceleration  vector  decoiiqwsed  into  an  axial  component, 
h^,  and  a  transverse  component,  ijj,  becomes 


b  =  [-<«)^(ro+  cos4))  cos<j)+  ( -gr)  sin<|)]  + 

[-0)2  (r' +  (z^-z)  cos<l))  sirwl)- (-g^)  cos<j)]^  +  (ro+ (z^-z)  cos<t)) 


do) 

dt 


211/J 


(1) 


Here,  the  transverse  component,  bj^,  is  a  magnitude  with  no  reference  to  direction  due  to 
the  axisymmetric  nature  of  the  capillary  tube.  This  form  of  the  acceleration  vector  does 
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not  include  the  Coriolis  acceleration  (see  Appendix  A)  which  would  be  embodied  in  the 
transverse  acceleration  component.  The  Coriolis  acceleration  is  induced  from  the  motion 
of  the  fluid  column  or  meniscus  in  the  axial  direction  in  the  capillary  tube.  Therefore 
prior  knowledge  of  the  meniscus  dynamics,  specifically  the  velocity,  is  required  to 
evaluate  the  Coriolis  acceleration.  The  one-dimensional  equation  of  motion  for  the 
meniscus  formulated  for  this  investigation  assumes  that  the  transverse  acceleration 
component  has  a  negligible  effect  on  the  dynamic  behavior  of  the  meniscus  and  therefore 
does  not  account  for  tangential  and  Coriolis  acceleration  effects. 

Equation  of  Motion 

The  simple  analytical  model  described  in  chapter  I  was  used  to  predict  the  motion 
of  the  meniscus  subjected  to  transient  accelerations.  The  analytical  formulation  of  the 
equation  of  motion  was  simplified  assuming  bulk  flow  only  in  the  axial  direction  of  the 
capillary  tube.  Effectively,  the  liquid  column  in  the  capillary  tube  was  assumed  to 
undergo  a  slug  flow. 

The  surface  force,  ,  was  formulated  by  accounting  for  the  cumulative  effect 

of  solid-liquid  intermolecular  forces  in  the  near  contact  line  region  and  the  wall  shear 
stress  associated  with  the  piedominandy  Poiseuille  flow  in  the  liquid  column.  This 
approach  has  been  used  since  the  work  of  West^®  and  Washbum^^  and  as  described  by 
Kafka  and  Dussan,  has  been  surprisingly  good  in  modelling  steady  capillary  flows 
when  Ca  and  Re  are  small.  The  contact  line  force  is  defined  assuming  a  spherical 
meniscus  and  a  constant  equivalent  or  apparent  contact  angle  such  that 
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=  a  (27ia)  cos6g  +  T  ^(27ta)  T] . 


(2) 


Note  that  such  an  assumption  implies  a  non-changing  apparent  contact  angle  under 
dynamic  conditions  (0  «  Q  ^). 

The  resulting  dimensional  form  of  the  momentum  equation  is 


2ocos6 

+ 

_8]Ln^ 

pa 

[  dt^  \  dt)  J 

pa" 

0)  ( t)  ^  cos(()/  ( -  -3- )  cos<j)  +  To)  +  gsincj) 


(3) 


The  transient  angular  velocity  is  of  the  form 


0)  ( t)  =  Bt, 

=  B{  tf-  t)  , 


tf 

2 


2 


(4) 


where 


U 
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Using  the  following  dimensionless  parameters, 


c  =  0  ^  c  ^  1 

0  =  0  i  1, 


(6) 


with  the  following  initial  conditions 

c  =  1,  -§  =  0,0  =  0, 


the  non-dimensional  form  of  the  momentum  equation  in  terms  of  a  transient 
dimensionless  meniscus  position  becomes 


-^sind) 

9pk 


a 

■  ReCa' 

+ 

16Cal 

\r  c?Cl 

l^oj 

2Bo 

de*  1  <*/  ; 

Bo 

Q  (0)  ^  cos<j) 


2 


jcOS<j)+  To 


+ 


-^sin4) 

9pk 


(7) 


where  the  angular  velocity,  Q ,  is  defined  by 


£3(0) 


O^0i  — 
2 

—  ^0^1. 
2 


(8) 


The  first  and  second  terms  of  eq.  (7)  represent  the  inertial  and  viscous  effects, 
respectively,  with  the  coefficients  formed  as  the  products  of  the  Reynolds  number, 
capillary  number.  Bond  number,  and  aspect  ratio.  The  Reynolds  and  capillary  numbers 
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are  referenced  to  a  characteristic  velocity,  which  is  the  maximum  attainable  velocity 
over  one-half  of  the  acceleration  cycle  period  assuming  t]  =  0  at  The  Bond 

number  is  referenced  to  the  peak  radial  acceleration  component, 

The  momentum  equation  can  be  easily  modified  to  include  a  constant  receding  or 
advancing  contact  angle,  0  differing  from  that  of  the  static  contact  angle,  0  resulting 
in 


g  COS0^j 
^  COS0^ 


sin4> 


a 

'  Re Ca' 

.^0. 

2Bo 

+ 

16Cal 

k-^1 

\d&j 

Bo 

[  d9j 

Q(0)^  cos<I)f[z^- 

cos<|)  +  rl  +  -2-  sind) 

VV  2  J 

/  ^pk 

(9) 


It  may  be  noted  that  the  general  form  of  the  dynamic  contact  angle  is  known  to  be 
functionally  dependent  on  the  contact  line  receding  and  advancing  velocity, 
0d  =  0<j(dC/ <d0).  Equation  (9)  can  be  modified  to  include  any  specific  form  for  the 
dynamic  contact  angle  and  solved  numerically. 

Results 

Equation  (9)  was  solved  for  the  dimensionless  meniscus  position,  C,  using  a 
central-difference  scheme.  Solutions  were  obtained  by  first  specifying  a  static  wicking 
height,  h^,  based  on  the  experimentally  determined  equivalent  contact  angle, 
0  5  =  16.7°.  Numerical  solutions  were  consistent  with  the  experimental  parameters  for 
Bo  =  0.3675,  1.47  and3.47xl0~®  ^  Ca.  ^  4 . 0  0  x  1 0  with  single  or  multiple 
acceleration  transient  cycles.  Numerical  solutions  were  also  calculated  for  various 
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advancing  contact  angles  in  order  to  address  deviations  between  e]q)erimental  results  and 
analytical  calculations  based  upon  the  assumption  that  the  contact  angle  is  equivalent  to 
the  static  contact  angle,  0  All  numerical  solutions  assumed  a  constant  recession  contact 
angle  equivalent  to  the  static  contact  angle,  0^  =  0^  =  16.7°. 

Experiments  were  performed  and  the  results  compared  with  solutions  to  the 
analytical  model.  Uncertainty  in  the  experimentally  obtained  dimensionless  meniscus 
position,  C,  and  time,  0,  was  on  the  order  of  C,  0  ±  0 . 02  based  on  the  accuracy  and 
resolution  of  the  experimentally  measured  quantities. 

Figure  2.6  compares  the  experimentally  obtained  dimensionless  meniscus  position 
to  the  analytical  results  at  varying  angles  of  inclination.  Also  shown  are  the  axial  and 
transverse  acceleration  components.  Experimental  parameters  in  this  figure  ranged  from 
3.47xl0-«  ^  Ca  ^  6.00xl0-«,  2.96x10-2  ^  Re  ^  5.13x10-2,  and 
5.71  X 10-2  ^  a/hg,  ^  9 . 89  X 10-2  ^  jjjg  inclination  angle  varied  from  30“  to  60“. 

The  experimentally  obtained  dimensionless  meniscus  position  agreed  weU  with  the 
analytical  calculations.  There  was  a  slight  deviation  from  the  predicted  values  due  to  the 
transverse  acceleration  component  as  shown  in  Fig.  2.6d  as  was  observed  in  chapter  one. 
A  small  Bond  number  tended  to  decrease  the  sensitivity  of  the  meniscus  position  to  the 
transverse  acceleration  component  as  a  result  of  the  small  C£q)illary  tube  diameter. 

Figure  2.7  compares  the  experimental  results  and  analytical  calculations  for 
4.00x10  2  ^  ca  ^  4.00x10“^  with  Bo  =  1.47  at  25“  inclination  and 

Bo  =  0.367  5  at  30“  inclination.  Table  2.1  shows  a  summary  of  the  dimensiotiless 
parameters  corresponding  to  Fig.  2.7. 
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Table  2.1 

Values  of  dimensionless  numbers. 


Bo 

Ca 

Re 

We 

alK 

<1) 

0.3675 

4  .  OOxlO'"* 

3.42 

1.37x10-5 

5 .71x10-5 

30° 

2.00x10-^ 

1.71 

3 .42x10-* 

1. 00x10-* 

0.855 

8.55x10-5 

4 .00x10-5 

0.342 

1.37x10-5 

6 .00x10-5* 

5.13x10-2 

3.08x10-’ 

1.47 

2.37x10-* 

4.05 

9 .60x10-* 

1.93x10-2 

25° 

1.18x10-* 

2.02 

2 .38x10-* 

5.92x10-5 

1.01 

5.98x10-5 

3.55x10-®** 

6.07x10-2 

2.15x10-’ 

*Fig.  2.6a 
“Fig.  1.6a 


As  shown  in  Fig.  2.7a,  there  was  no  significant  deviation  between  the  predicted 
and  observed  meniscus  position  for  Bo  =  1.47.  For  Bo  =  0.3675, as  shown  in  Fig. 
2.7b,  the  predicted  and  observed  meniscus  position  agreed  well  for  meniscus  recession 
with  an  increase  in  the  retardation  in  meniscus  position  corresponding  to  an  increase  m 
the  capillary  number.  For  the  advancing  meniscus,  there  was  a  retardation  in  the 
observed  meniscus  position,  in  response  to  the  acceleration  transient,  from  that  of  the 
predicted  values.  The  retarded  position  for  the  advancing  meniscus  was  not  evident  for 
Ca  =  4.00  xlO"^  but  became  significant  for  Ca  ^  1 . 00  xlO'^  and  increased  with 
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increasing  capillary  number. 

The  variation  between  the  predicted  and  observed  advancing  meniscus  velocity  is 
in  part  due  to  the  dependence  of  the  contact  angle  on  capillary  number  and/or  Reynolds 
number.  Figure  2.8  compares  the  anal5^ical  calculations  for  a  16.7°  receding  contact 
angle  while  varying  the  advancing  contact  angle.  Also  shown  in  Fig.  2.8  are  the  transient 
acceleration  fields  decomposed  into  the  axial  and  transverse  components  at  the  pivot  point 
of  the  capillary  tube.  The  experimental  results  agree  very  well  with  that  of  the  analytical 
calculations  when  a  change  in  the  advancing  contact  angle  is  taken  into  account. 

Hoffinan*^  and  Phillips  and  Riddiford*'*  discuss  systems  in  which  the  inteifacial 
forces  between  the  solid  and  the  liquid  change  when  flow  occurs  under  conditions  where 
inertia  is  assumed  to  be  negligible.  For  very  low  velocities  the  £q)parent  contact  angle  is 
usually  independent  of  velocity,  but  as  the  inteifacial  velocity  is  increased  above  a  critical 
velocity,  the  apparent  contact  angle  will  change  until  a  limiting  value  is  reached.  Phillips 
and  Riddiford  postulated  that  there  is  a  change  in  the  interfacial  forces  as  the  critical 
velocity  is  reached  thereby  resulting  in  a  shift  in  the  qjparent  contact  angle.  Figures  2.8b 
and  8c  show  that  the  advancing  velocity  agrees  well  with  the  analytical  calculations  up 
to  0  -  0  .75,  corresponding  to  an  inflection  in  the  transverse  acceleration  component. 
The  meniscus  velocity  at  0-0.75  suddenly  decreases  (from  -3  mm/s  for 

Ca  =  2.00 xlO”^  and  -6  mm/s  for  ca  =  4.00 xlO”'®)  and  subsequently  increases 
again  resulting  in  a  retarded  meniscus  position  compared  to  the  analytical  solution.  The 
sudden  decrease  in  the  meniscus  velocity  at  time  0  -  0.75  and  subsequent  increase  may 
be  due  to  the  combined  effects  of  the  transverse  acceleration  component  and  the  Coriolis 
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acceleration  which  is  dependent  on  the  velocity  of  the  meniscus.  In  addition,  the 
analytical  formulation  takes  into  account  the  opposing  viscous  force  at  the  tube  wall  and 
general  inertial  effects  but  does  not  take  into  account  the  opposing  viscous  force  at  the 
contact  line  region  and  a  possible  increase  in  the  radius  of  curvature  of  the  meniscus  as 
a  result  of  inertial  effects.  Analytical  calculations  and  experimental  results  show  that  the 
motion  of  the  meniscus  along  the  contact  line  and  fluid  along  the  tube  wall  is  retarded 
due  to  viscous  effects.  However,  the  motion  of  the  central  core  of  fluid,  due  to  inertia, 
likely  results  in  the  "pushing  forward"  of  the  central  portion  of  the  meniscus  relative  to 
the  contact  line  yielding  an  increased  radius  of  curvature  and  therefore  an  increased 
contact  angle.  This  change  in  contact  angle  further  retards  the  meniscus  motion. 

In  addition,  with  the  removal  of  the  acceleration  transient,  there  is  a  rt^id  but 
noticeable  relaxation  time  ( «*4-5  seconds)  when  the  contact  angle  returns  to  the  static 
value  and  the  column  of  alcohol  returns  to  the  initial  static  wicking  height.  The  existence 
of  a  relaxation  time  agrees  with  observations  by  Phillips  and  Riddiford.*^ 

Figure  2.9  compares  the  analytically  calculated  advancing  contact  angle  required 
to  match  experimental  results  to  the  results  presented  by  Hoffman’^.  The  shifting  function 
described  by  Hof&nan  was  determined  to  be  ^jproximately  F  (6  g)  =2.5x10“^.  Using 
a  capillary  number  as  defined  in  this  investigation  (really  an  average  capillary  number), 
Ca  +  F(0  g)  was  plotted  against  the  advancing  contact  angle.  The  deviation  from  that 
of  Hoffinan’s  results  may  be  in  part  due  to  the  fact  that  Hoffman  is  characterizing  the 
shape  of  the  liquid-air  interface  for  viscous  fluids  (i.e.  oils)  in  motion  when  only  the 
viscous  and  interfacial  forces  are  important.  This  investigation  considers  motion  of  the 
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liquid-air  interface  for  ethyl  alcohol  as  a  result  of  acceleration  induced  forces  which 
results  in  inertial  as  well  as  the  viscous  and  inteifacial  forces. 

Figure  2.10  shows  the  percent  deviation  in  meniscus  position  from  the  inviscid 
solution,  at  the  end  of  a  cycle,  relative  to  the  viscous  term,  Ca/Bo.  The  effect  of 
increasing  the  angle  of  inclination  resulted  in  a  reduction  in  the  percent  deviation  from 
the  inviscid  solution.  As  shown  in  Fig.  2.10,  the  experimental  results  showed  an 
increased  deviation  over  the  anal)rtical  solutions  which  appears  as  a  vertical  shift  from  the 
analytical  results.  This  shift  between  the  experimental  results  and  analytical  predictions 
provides  credibility  to  the  hypothesis  that  inertial  forces  may  result  in  a  change  in  the 
contact  angle  and  that  viscous  forces  alone  are  not  sufficient  to  induce  the  observed 
increase  in  the  deviation  (or  attenuation)  from  the  inviscid  prediction. 

Figure  2.11  compares  analytical  calculations  and  experimental  results  for  three 
cycles  m  the  transient  acceleration  for  So  =  0.367  5,  Ca  =  2.00x10"^  and 
4.00x10"^,  =  2 .  ig,  4 . 2g,  8 . 4g,  and  at  a  30°  inclination  angle.  As  shown 

in  Fig.  2.11a,  there  is  an  attenuation  and  phase  shift  in  meniscus  position  due  to  viscous 
effects  from  an  inviscid  or  low  capillary  number  solution  as  the  capillary  number  is 
increased.  Again,  there  is  good  agreement  between  the  observed  and  predicted  receding 
meniscus  position.  The  deviation  between  the  observed  and  predicted  advancing  meniscus 
position  is  apparent,  however,  when  an  advancing  contact  angle  of  30°  is  accounted  for 
in  the  numerical  solution.  The  predicted  and  observed  meniscus  positions  agree  very 
weU.  As  shown  in  Fig.  2.11b,  there  is  an  increased  attenuation  from  the  initial  meniscus 
height  and  a  slight  phase  shift  with  increasing  peak  radial  acceleration. 
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Fig.  2.1  Schematic  of  centrifuge. 
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Fig.  2.2  Plot  of  the  centrifuge  angular  velocity  test  cycle  (f=  0.0015  Hz). 
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acceleration, 


a) 


Fig.  2.3  Experimentally  generated  acceleration  field  (accelerometer  output): 
a)  r  =  0 . 0015  Hz,  single  cycle, 
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Fig.  2.3(contd)  b)  jf  =  0.1  Hz,  multiple  cycle, 
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Fig.  2.3(contd)  c)  f  =  0.05  Hz,  multiple  cycle, 
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Fig.  2.3(contd)  d)  =  0 . 025  Hz,  multiple  cycle, 
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Fig.  2.3(contd) 
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e)  f  =  0.01  Hz,  multiple  cycle. 
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Fig.  2.4  Experimental  apparatus:  a)  sealed  test  cell 
and  b)  schematic  of  experimental  apparatus  as  mounted  on  the  centrifuge. 
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a) 


Fig.  2.5  Capillary  tube  orientation:  a)  pivot  point  location,  directional  coordinates  and 
b)  coordinate  systems  as  referenced  to  the  centrifuge  accelerometer  and  capillary  tube. 
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Fig.  2.6  Experimental  data  compared  to  analytical  calculations 
for  ethyl  alcohol  at  So  =  0.3675: 

a)  d)  =  30®,  ca  =  6.00x10-®,  Re  =  5.13x10-2,  a/h„  =  5.71x10 


Fig.  2.6(contd) 

b)  4)  =  40°,  Ca  =  4.67  xlO'®,  Re  =  3.99x10-2  ^  7 .34x10 


acceleration, 


acceleration, 


Fig.  2.6(contd) 

d)  <1)  =  60°,  Ca  =  3.47x10-®,  Re  =  2.96x10-2,  a/h^  =  9.89x10 
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Ca  =  2.37  X  1  O”'*,  numerical 
Ca  =  1 .1  8  X  1  O"**,  numerical 
Ca  =  5.92  X  1  0“^,  numerical 
Ca  =  2.37  X  1  experimental 
Ca  =  1  .1  8  X  1  0“'*,  experimental 
Ca  =  5.92  X  1  0“^,  experimental 
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2.7  Experimental  data  compared  to  analytical  calculations 
for  varying  cjqjillary  number 
a)  d)  =  25°,  =  1.47, 


acceleration, 


-  9d  =  1  6.7°  advancing,  numerical 

- 9d  =  24°  advancing,  numerical 

o  ■  experimental  results 
-S--  bz 
-0--  bR 


Fig.  2.8(contd)  b)  Ca  =  2.00x10"*, 


acceleration 


0d  advanc ing,  deg. 


Fig.  2.9  Conyarison  of  Hoffinan’s*^  results  and  advancing  contact  angle 
required  to  match  experimental  results. 
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Fig.  2.10  Percent  deviation  from  the  inviscid  solution  at  the  end  of  a  cycle. 
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in  vise  id  soil.,  0d  =  1  6.7°  advanc  ing 

Co  =  4.00  X  1  O”"*,  9d  =  1  6.7°  advancing,  numerica 

Ca  =  2.00  X  1  O""*,  9d  =  1  6.7°  advancing,  numerica 

Co  =  4.00  X  1  0"'*,  experimental 

Ca  =  4.00  X  1  0"'*,  9d  =  30°  advancing,  numerical 


Fig.  2.11  Experimental  data  compared  to  analytical  calculations  for 
multiple  acceleration  cycles  =  30°,  Bo  =  0.3675): 
a)  varying  capillary  number,  -  A  .2 g, 


CHAPTER  m 


DYNAMIC  BEHAVIOR  OF  A  HEATED  MENISCUS 
SUBJECTED  TO  A  TRANSIENT  ACCELERATION  HELD 

Abstract 

This  study  e^erimentally  and  analytically  investigates  the  effects  of  acceleration 
transients  on  the  dynamic  response  of  a  heated  meniscus  in  a  capillary  tube.  A  capillary 
tube  partially  filled  with  ethyl  alcohol  was  moimted  on  a  centrifuge  to  observe  the 
dynamic  response  of  a  heated  meniscus  subject  to  an  acceleration  transient.  An  initial 
recession  of  the  meniscus  due  to  heat  iiq^ut  was  treated  as  a  change  in  contact  angle  due 
to  thermocapillary  effects  as  a  result  of  non-isothermal  heating.  Experiments  and 
analytical  calculations  were  for  a  Bond  mrmber  of  1.47  and  a  capillary  number  of 
3.00x10"®.  Experimental  data  of  the  meniscus  recession  and  advance  agreed  well  with 
a  one-dimensional  equation  of  motion  accounting  for  the  change  in  contact  angle  due  to 
thermoc^illary  effects. 


Experimental 
Acceleration  Field 

Transient  acceleration  fields  were  generated  using  a  2.4-m  diameter  centrifuge 
table  (Fig.  3.1)  rotating  with  a  time  variant  angular  velocity.  Angular  velocity  transients 
were  generated  at  a  cyclic  frequency,  f  =  0.0015Hz,  using  a  signal  or  waveform 
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generator  as  a  control  input  to  the  centrifuge.  The  angular  velocity  was  increased  linearly 
to  a  peak  value  at  the  midpoint  of  the  cycle  and  subsequently  decreased  for  two  cycles. 
A  tri-axis  accelerometer  (Columbia  Research  Laboratories,  Inc.)  was  used  to  monitor  the 
time  variant  acceleration  components  in  a  cartesian  reference  frame  affixed  to  the 
centrifuge  table.  For  this  investigation,  the  peak  angular  velocity  was  specified  such  that 
a  4.2-g  ±  0.1-g  peak  radial  acceleration  was  generated  =  41.20  m/s^). 

Capillary  Test  Cell 

A  sealed  test  cell  was  fabricated  of  delrin  to  contain  an  ethyl  alcohol  working  fluid 
and  a  1.0-mm  diameter  capillary  tube  as  shown  in  Fig.  3.2a.  The  capillary  tube  was  held 
in  place  using  spring  loaded  positioning  clamps.  A  transparent  window  allowed  viewing 
of  the  capillary  tube  during  experimental  testing.  A  1  cm  long  transparent  gold  film  was 
vapor  deposited  on  the  top  exterior  surface  of  the  capillary  tube  and  served  as  a  heater.^* 
A  Kepco  precision  dc  power  supply  (Kepco,  Inc.)  was  used  to  power  the  gold  film  heater. 

Exposed  junction  copper-constantan  type  thermocouples  were  attached  at  a  spacing 
of  approximately  1  cm  to  the  exterior  surface  of  the  capillary  tube  with  a  grqjhite  slurry 
to  maintain  a  thermal  bond  with  an  outside  coating  of  high  temperature  ceramic  to 
maintain  strength  integrity.  Additional  thermocouples  were  placed  in  the  liquid  reservoir 
and  vapor  space  (properties  of  the  ethyl  alcohol  were  assumed  to  be  that  of  the  bulk 
reservoir  fluid  corresponding  to  an  experimental  temperature  of  30“C).  Thermocouple 
signals  were  linearized  and  amplified  to  a  0  -  5  volt  signal  on  the  centrifuge  to  eliminate 
slip  ring  noise  and  error  due  to  dissimilar  metal  junctions.  The  data  acquisition  system 
(including  signal  conditioner,  slip  ring,  A/D  converter,  etc.)  and  thermocouples  were 
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calibrated  over  a  temperature  range  of  25  -  TO^C  using  a  high  precision  temperature  probe 
(Hart  Scientific  5690)  to  within  ±  l^C. 

The  test  cell  and  capillary  tube  were  cleaned  using  potassium  hydroxide  and 
thoroughly  rinsed  with  ethyl  alcohol.  The  test  cell  was  filled  and  subsequendy  sealed 
resulting  in  a  fill  that  had  a  combination  of  air,  fluid  vapor,  and  liquid.  The  sealed  test 
cell  containing  a  glass  capillary  tube  and  reservoir  was  positioned  on  the  centrifuge  as 
shown  in  Fig.  3.2b.  The  test  cell  was  moimted  to  a  motorized  optics  rotation  stage 
(Newport)  such  that  the  capillary  tube  was  allowed  to  pivot  about  its  center  of  mass.  The 
dynamic  response  of  the  meniscus  to  the  transient  acceleration  field  was  observed  using 
an  8-mm  format  CCD  camera  (Sony,  30  frames/s)  mounted  adjacent  to  the  test  cell. 
Video  images  were  used  to  provide  a  record  of  the  meniscus  height  relative  to  the 
reservoir  meniscus,  t|  ,  as  a  function  of  time  to  within  ±1.0  mm.  An  equivalent  static 
contact  angle,  0^,  (with  and  without  the  addition  of  heat)  was  calculated  from  the 
experimentally  determined  static  wicking  height,  for  conditions  when  the  transverse 
acceleration  component  was  zero  (h^pg  -  (2a /a)  cosB^). 

Analytical  Formulation 

The  goals  of  the  analytical  formulation  were  to  expand  the  mathematical 
formulation  described  in  chapters  I  and  n  to:  (i)  mathematically  describe  the  dynamic 
behavior  of  a  heated  meniscus  as  a  result  of  the  temporal  acceleration  induced  forces 
using  a  simplified  one-dimensional  model  accounting  for  changes  in  the  dynamic  contact 
angle  due  to  heat  addition;  (ii)  address  the  altered  dynamic  behavior  of  the  heated 
meniscus  by  accounting  for  thermocapillary  effects  through  a  change  in  the  dynamic 
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contact  angle;  and  (iii)  compare  experimental  results  and  analjrtical  calculations. 

Acceleration  Vector 

Acceleration  measurements  were  obtained  with  regard  to  a  three-dimensional 
cartesian  non-inertial  reference  frame  affixed  to  the  centrifuge  as  described  in  chapter  one. 

The  acceleration  field  was  then  referenced  to  fixed  locations  on  the  capillary  tube  and 
subsequently  decomposed  into  transverse  and  axial  acceleration  components  relative  to 
the  capillary  tube  as  shown  in  Fig.  3.3. 

The  resulting  transient  acceleration  vector  decomposed  into  an  axial  component, 
jbj.,  and  a  transverse  component,  becomes 

b  -  [-«^(ro+  (Zj-z)  cos<|))  cos<|)+  ( -gr)  sinrj)] ^2  + 

[-Ci)2(ro+  (z^-z)  cos<{))  sin<|)-  (-gr)cos4)]^  + 

(1) 

Here,  the  transverse  component,  Jd^j,  is  a  magnitude  with  no  reference  to  direction  due  to 
the  axisymmetric  nature  of  the  capillary  tube. 

Equation  of  Motion 

A  simple  analj^ical  model  taking  into  account  changes  in  the  dynamic  contact 
angle  upon  recession  and  advance  was  used  to  predict  the  motion  of  the  meniscus 
subjected  to  transient  accelerations.  The  analytical  formulation  of  the  equation  of  motion 
was  simplified  assuming  bulk  flow  only  in  the  axial  direction  of  the  capillary  tube. 
Effectively,  the  liquid  column  in  the  capillary  tube  was  assumed  to  undergo  a  slug  flow. 
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The  contact  line  force  is  defined  assuming  a  spherical  meniscus  and  a  constant  receding 
or  advancing  contact  angle,  d^,  differing  from  that  of  the  static  contact  angle,  0^. 

The  resulting  dimensional  form  of  the  momentum  equation  is 


+  T1 


2ocos6 
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_  dt^  [dtj 

dtJ 
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The  transient  angular  velocity  is  of  the  form 
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Using  the  following  dimensionless  parameters, 
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with  the  following  initial  conditions 


the  non-dimensional  fonn  of  the  momentum  equation  in  tenns  of  a  transient 
dimensionless  meniscus  position  becomes 
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where  the  angular  velocity,  Q ,  is  defined  by 
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The  first  and  second  terms  of  eq.  (6)  represent  the  inertial  and  viscous  effects, 
respectively,  with  the  coefficients  formed  as  the  products  of  the  Reynolds  number, 
capillary  number.  Bond  number,  and  aspect  ratio.  The  Reynolds  and  cJ5)illary  numbers 
are  referenced  to  a  characteristic  velocity,  which  is  the  maximum  attainable  velocity 
over  one-half  of  the  acceleration  cycle  period  assuming  q  =  0  at  The  Bond 

number  is  referenced  to  the  peak  radial  acceleration  component,  gp^. 


ThermocapiUary  Effects 

With  addition  of  heat  into  the  meniscus  region,  it  has  been  speculated  that  thetmocapillary 
stresses  cause  an  increase  in  the  apparent  contact  angle  as  discussed  by  Pratt  and 
HaUinan.^®  The  thermocapillary  force  defined  by  Pratt  and  Hallinan  can  be  cast  as  a 
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function  of  the  axial  temperature  gradient  along  the  capillary  tube  and  the  gradient  of 
surface  tension  to  temperature  such  that 


da  dT 


da 

dT 


(8) 


Ar(27ia) 


assuming  a  constant  waU  temperature  gradient.  Using  this  approximation,  it  follows  that 
a  change  in  static  contact  angle  as  a  result  of  heat  addition  can  be  accounted  for  where 


0(27ta)C0se^_  Aeat  = 


O(27ia)cos0^-||  U^^(27ra) 

a  {2na)  cos0  -^Ar(27ta) 
s  Sr 


4^Ar  -  o  (cos0  _-cos0„ 
dT 


heat* 


(9) 


If  it  is  further  assumed  that  heat  addition  produces  a  change  in  contact  angle  from  the 
static  value  to  a  constant  dynamic  contact  angle  such  that  0  ^  =  0^  then  the  derived 
analytical  form,  eq.  (6),  can  be  solved  to  determine  the  meniscus  dynamic  response  to  an 
acceleration  transient  with  the  addition  of  heat  provided  that  0^  ^  known.  This 

analysis  assumes  negligible  buoyancy  effects  as  discussed  by  Pratt  and  HaUinan.“ 

Results 

Equation  (6)  was  solved  for  the  dimensionless  meniscus  position,  using  a 
central-difference  scheme.  Numerical  solutions  were  consistent  with  the  experimental 
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parameters  for  Bo  =  1.47  and  Ca  =  3.00x10"®  with  two  acceleration  transient 
cycles.  For  the  case  with  no  heat  addition  and  due  to  the  slow  acceleration  transient, 
solutions  were  obtained  by  specifying  the  experimentally  inferred  static  contact  angle  (see 
Appendix  B)  to  be  equivalent  to  the  dynamic  contact  angle,  0^  =  0^^  =  16.7®.  With 

the  addition  of  heat  to  the  capUlary  tube  there  was  an  immediate  recession  of  the  contact 
line  to  a  new  static  wicking  height  corresponding  to  an  altered  contact  angle,  0^ 
accounting  for  heat  addition.  The  value  of  this  contact  angle  was  calculated  from  the 
experimental  results  as  was  the  static  contact  angle  for  the  unheated  case  (also  see 
Appendix  B).  Numerical  solutions  for  the  case  with  heat  addition  assumed  the  static 
contact  angle,  0  ^  =  16.7®,  and  a  constant  dynamic  contact  angle  equivalent  to  the 
e^)erimentally  inferred  contact  angle  at  the  retarded  wicking  height  due  heat  addition 
such  that  0j  =  0s^^aat  =  35.6°.  Fluid  properties  ( o ,  p ,  p )  were  referenced  to  the 
bulk  fluid  in  the  reservoir,  =  30  °C.  The  temperature  gradient,  AT,  along  the 
capillary  tube  was  determined  using  thermocouple  data  and  was  assumed  to  be 
representative  of  the  A  r  between  thermocouples  tc2  and  tc3  just  below  the  heater. 

Ejqwriments  were  performed  at  a  fixed  angle  of  inclination  of  30°  while  applying 
heat  to  the  top  exterior  of  the  capillary  tube.  The  acceleration  transient  was  not  applied 
until  thermal  equilibrium  was  obtained  within  the  test  cell.  Three  levels  of  heat  inputs, 
including  line  losses,  were  1.0 , 3.5  W,  and  4.0  W  ±  0.5  W.  The  dynamic  response  of 
the  heated  meniscus  was  compared  with  solutions  to  the  analytical  model.  Uncertainty 
in  the  experimentally  obtained  dimensionless  meniscus  position,  Q,  and  time,  0,  was  on 
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the  order  of  C ,  0  ±  0.02  based  on  the  accuracy  and  resolution  of  the  experimentally 
measured  quantities. 

Figure  3.4  shows  the  temperature  profile  along  the  capillary  tube  and  within  the 
test  cell  for  each  of  the  three  heat  inputs.  Also  shown  are  the  transient  axial,  and 
transverse,  jb^,  components  of  the  acceleration  vector  at  the  pivot  point  of  the  capillary 
tube  for  two  transient  cycles.  For  clarification  in  Fig.  3.4c,  the  arrows  at  the  uppermost 
part  of  the  graph  show  the  dimensionless  time  for  which  the  dimensionless  meniscus 
position  corresponds  to  the  indicated  thermocouple  location.  The  temperature  variations 
during  the  experimental  cycles  are  due  to  movement  of  the  meniscus  in  the  capillary  tube 
as  well  as  the  orientation  of  the  reservoir  surface  with  respect  to  the  capillary  tube.  As 
the  liquid  column  moves  either  out  of  or  into  the  capillary  tube,  there  is  a  corresponding 
increase  or  decrease  in  the  rate  of  energy  storage  within  the  capillary  tube.  This  change 
in  energy  storage  results  in  a  change  in  the  temperature  profile  of  the  capillary  tube. 
Also,  as  the  acceleration  vector  varies  in  magnitude  and  direction,  the  reservoir  surface 
remains  in  an  orientation  normal  to  the  acceleration  vector  and  therefore  is  continually 
in  motion  as  is  the  liquid  column.  The  location  of  the  reservoir  surface  relative  to  the 
proximity  of  the  thermocouples  will  influence  the  local  capillary  mbe  temperature. 

The  variation  in  the  temperature  profile  of  the  capillary  tube  as  a  result  of  the 
motion  of  both  the  liquid  column  and  reservoir  is  depicted  in  Fig.  3.5.  Figure  3.5  shows 
photographs  of  the  meniscus  and  reservoir  position  and  the  temperature  history  for 
selected  times  during  the  first  cycle  with  a  heat  input  of  4.0  W,  Fig.  3.4c.  The 
photographs  were  taken  at  times  corresponding  to  ©  =  0.0735,  0.1605,  0.2535,  0.4320, 
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0.7065, 0.8280, 0.9090  or  to  the  maximum  non-dimensional  time  indicated  on  the  0  axis 
of  the  acconq)anying  temperature  graph.  From  Fig.  3.5b,  as  the  meniscus  passes  between 
thermocouples  tc2  (located  just  below  the  heater)  and  tc3  there  is  an  increases  in 
temperature  in  tc2  and  tc3  and  a  decrease  in  the  reservoir  temperature.  As  the  meniscus 
passes  between  thermocouples  tc3  and  tc4,  Fig.  3.5c,  the  reservoir  temperature  stabilizes 
with  continued  increases  in  temperature  of  tc2  and  tc3.  It  may  be  noted  that 
thermocouple  tc4,  located  in  close  proximity  of  the  reservoir  surface  at  the  start  of  the 
cycle,  increases  slightly  in  tenqjerature  and  stabilizes  as  the  reservoir  is  reoriented  below 
the  thermocouple.  There  is  also  an  inflection  in  the  transient  temperature  response  of  tc2 
and  tc3  as  the  transverse  acceleration  component  passes  through  a  minimum.  As  the 
acceleration  reaches  a  maximum  and  begins  to  decline.  Fig  3.5d  and  3.5e,  the  rate  of 
change  of  tc2  and  tc3  decreases  toward  a  peak  temperature  and  there  is  a  sudden  increase 
in  vapor  temperature.  At  this  point  in  time  the  meniscus  position  is  at  a  minimum  and 
conduction  of  heat  down  the  c^illary  tube  is  ^proaching  steady  state  conditions  with  a 
reduction  in  the  energy  storage  (due  in  part  to  the  absence  of  the  column  of  liquid  in  the 
capillary  tube)  in  the  capillary  tube  itself.  As  a  result,  the  excess  energy  is  going  into  the 
v^Ktr  causing  a  temperature  rise.  As  the  acceleration  is  decreased  for  the  remainder  of 
the  cycle,  the  return  path  of  the  meniscus  and  reservoir  surface  motion.  Figs.  3.5  e-g, 
result  in  the  reversal  in  the  transient  temperature  response. 

Figure  3.6  compares  experimental  results  and  analytical  calculations  for  each  of 
the  three  heat  inputs  for  two  acceleration  transient  cycles.  Also  shown  are  the 
components  of  tiie  acceleration  vector.  There  was  no  significant  difference  in  the 


89 


measured  dimensionless  meniscus  position  due  to  the  heat  inputs.  Analytical  solutions 
consisted  of:  (i)  maintaining  a  constant  contact  angle  such  that  Q  =  0  ^  =  16 . 7  (ii) 
assuming  a  change  in  contact  angle  due  the  thermocapillary  effects  such  that 
=  35.6°,  0g  =  16.7°;  (iii)  assuming  a  change  in  contact  angle  due  the 

theimocapillary  effects  such  that  0^;  =  40.0°,  0^  =  16.7°;  and  (iv)  the  dynamic 
contact  angle,  =  0^  was  constant  for  the  three  heat  inputs.  Analytical  solutions 
consistent  with  the  assun^tion  that  there  is  a  change  in  contact  angle  due  to 
thermocapillary  effects  showed  good  correlation  with  the  experimental  results  for  the  three 
heat  inputs  for  the  case  where  0^=  35.6°.  Solutions  for  0  j  =  40 . 0  °  resulted  in  a 
retarded  meniscus  height  compared  to  the  experimental  results. 

The  results  discussed  so  far  have  assumed  that  thermocapiUary  effects  due  to 
heating  of  the  meniscus  will  result  in  a  change  in  contact  angle  and  that  this  contact  angle 
can  subsequently  be  used  to  predict  the  meniscus  motion  with  heating.  Equation  (9)  can 
be  verified  using  experimental  data  to  evaluate  AT,  0^, ,  and  =  0^  and 

standard  property  tables  for  o  and  dc/dT.  An  average  AT  was  calculated  from 
experimental  data  using  the  difference  between  the  thermocouples  tc2  and  tc3.  Table  3.1 
summarizes,  the  information  required  to  evaluate  the  terms  in  eq.  (9)  for  the  heat  inputs 
of  1.0,  3.5,  and  4.0  W.  For  a  heat  input  of  4.0  W,  the  error  between  the  two  terms  was 
within  5%  while  for  a  heat  input  of  1.0  W,  the  error  was  76%.  This  discrepancy  may  be 
in  part  due  to  the  assumption  that  the  fluid  properties  were  evaluated  at  the  reservoir 
temperature  as  well  as  error  in  the  experimentally  inferred  contact  angle,  0^  In 
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addition,  the  defined  ten^rature  gradient,  AT,  is  a  critical  value  for  which 
thermocapillary  stresses  will  result  in  an  altered  contact  angle.  In  general,  an  increase  in 
temperature  will  decrease  the  density  and  surface  tension.  The  combined  effects  of  a 
decrease  in  density  and  a  reduced  surface  tension  for  a  fixed  contact  angle,  assuming 
negligible  thermocapillary  stresses,  may  result  in  a  recession  in  the  static  wicking  height 
due  to  heat  input.  If  in  fact  density,  surface  tension,  and  contact  angle  change  due  to  the 
combined  effects  of  an  increase  in  temperature  and  thermocapillary  stresses  (or  a  critical 
A  T)  as  a  result  of  an  increased  heat  input,  there  may  be  situations  for  which  the  static 
wicking  height  may  not  be  altered  firom  the  case  where  thermocapillaiy  effects  were 
assumed  to  be  negdigible.  This  may  be  why  the  combined  assumptions  of  a  constant 
0g  and  evaluation  of  fluid  properties  at  SOT  resulted  in  a  good  correlation  between 

analytical  calculations  and  experimental  results. 
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Table  3.1 

Summary  of  thermoce^illary  terms  in  eq.  (9). 


o,  (@30°C),  (N/m) 

0.02189 

-1^,  [2Q°  ^  T  ^  100°) 
dT 

-9 .15x10"^ 

(N/m-°C) 

0  s 

16.7 

0d  =  0s.  ftsac  (deg.) 

35.6 

AT  CO 

-33.3  ±  3.3 

-28.9  ±  2.7 

-7.9  ±  0.5 

AT  (N/m) 

3 . 02x10'^ 

±  0.27x10-3 

2 .64x10-3 
±  0.25x10-3 

7 .2  xlO-^ 

±  0.5x10-“* 

o(cOSe^-COS0^,  i:sat) 

3.17x10-3 

(N/m) 

0(W) 

4.0 

in 

(n 

o 

H 
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Fig.  3.2  Experimental  apparatus:  a)  sealed  test  cell 
and  b)  schematic  of  experimental  apparatus  as  mounted  on  the  centrifuge. 
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a) 


Fig.  3.3  Capillary  tube  orientation:  a)  pivot  point  location,  directional  coordinates  and 
b)  coordinate  systems  as  referenced  to  the  centrifuge  accelerometer  and  capillary  tube. 
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CHAPTER  IV 


SUMMARY  AND  CONCLUSIONS 

In  summary,  the  dynamic  response  of  an  unheated  and  heated  meniscus  within  a 
capillary  tube  subjected  to  single  and  multiple  cycle  transient  acceleration  induced  forces 
was  studied  experimentally  and  analytically.  A  capillary  tube,  contained  within  a  test  cell, 
was  partially  Med  with  either  water  or  ethyl  alcohol  and  mounted  on  a  centrifuge  for 
visual  observation  of  the  dynamic  response  of  the  meniscus  subjeaed  to  the  transient 
acceleration  field.  Experimental  parameters  were  specified  such  that  the  capillary  and 
Bond  numbers  were  of  sufficient  magnitude  to  account  for  viscous,  inertial,  and  combined 
axial  and  transverse  acceleration  effects  on  the  unheated  and  heated  meniscus.  A  one¬ 
dimensional  mathematical  model,  assuming  a  spherically  sht^jed  meniscus,  was 
formulated  to  predict  the  dynamic  response  of  the  meniscus  to  an  acceleration  transient 
and  compared  with  experimental  results.  For  the  case  in  which  the  meniscus  was  heated, 
thennoc£^illary  effects  were  modeled  by  accounting  a  change  in  contact  angle  due  to 
heating  and  incorporated  in  the  one-dimensional  equation  of  motion. 

Comparison  of  the  experimental  results  and  analytical  calculations  for  the  case  of 
an  unheated  meniscus  showed  that  a  combined  axial  and  elevated  transverse  acceleration 
significantly  altered  the  dynamic  response  of  the  meniscus  and  resulting  in  an  increased 
meniscus  velocity  and  subsequently  lower  meniscus  position.  Systems  characterized  by 
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large  Bond  numbers,  greater  than  1.5,  as  distinguished  by  the  larger  capillaiy  tube 
diameters  were  particularly  sensitive  to  an  elevated  transverse  acceleration  component 
greater  than  2g.  It  is  speculated  that  this  dynamic  response  in  meniscus  position  is  due 
to  a  reduced  overall  contact  line  force  in  the  near  contact  line  region  resulting  from  the 
meniscus  being  moved  in  a  transverse  direction  and  distended  in  the  axial  direction. 
Systems  typified  by  a  low  Bond  number,  less  than  1.5,  tended  to  be  insensitive  to  a 
transverse  acceleration  component  due  to  the  small  capillary  tube  diameter. 

For  systems  typified  by  the  combination  of  a  low  Bond  number,  less  than  1.5,  and 
high  capillary  number,  greater  than  10 there  was  a  significant  change  in  the  dynamic 
response  of  the  meniscus.  There  was  a  significant  amoimt  of  retardation  in  the  meniscus 
position,  when  compared  to  an  inviscid  solution,  due  to  a  combination  of  viscous  forces, 
acceleration  induced  inertial  forces,  and  interfacial  forces.  It  may  be  noted  that  variations 
in  either  the  receding  contact  angle  or  the  advancing  contact  angle  are  indicative  of 
changes  in  any  or  all  of  these  forces  described. 

Hysteresis  and  sticking  of  the  advancing  meniscus  exhibited  with  water  may  be 
attributed  to  such  phenomena  as:  (i)  surface  roughness;  (ii)  contamination;  or  (iii)  draining 
or  evaporating  of  the  thin  film  left  ftom  the  receding  contact  line  resulting  in  the 
meniscus  advancing  onto  a  dry  surface. 

When  the  meniscus  was  heated,  without  an  induced  acceleration,  there  was  an 
initial  recession  of  the  meniscus  that  was  accounted  through  thermocapillary  stresses 
bringing  about  a  change  in  the  contact  angle.  Experimental  results  and  analytical 
calculations  for  a  Bond  number  of  1.47  and  a  capillary  number  of  3 . 00  xlO'®  agreed 


108 


well  when  accounting  for  the  change  in  contact  angle,  due  to  theimocapillaiy  effects,  in 
the  dynamic  contact  angle  term  in  the  equation  of  motion. 

Experiments  and  analytical  calculations  were  not  performed  for  cases  in  which 
viscous,  inertial,  and  combined  axial  and  transverse  acceleration  effects  aU  contribute  to 
the  dynamic  response  of  the  heated  meniscus.  The  combined  effects  of  heating,  viscous, 
inertial,  and  combined  axial  and  transverse  accelerations  may  result  in  a  meniscus 
dynamic  response  differing  significantly  from  observations  noted  in  this  investigation. 
In  addition,  the  extent  that  the  near  contact  line  region  affects  the  transient  behavior  of 
the  meniscus  was  not  addressed  in  this  investigation  except  through  changes  in  the 
dynamic  contact  angle. 

Some  suggested  areas  of  future  research  as  a  direct  result  of  this  investigation 
include  the:  (i)  meniscus  dynamic  response  to  steady  state  heating  with  a  combination  of 
low  Bond  number  and  high  capillary  number;  (ii)  meniscus  dynamic  response  to  increased 
heat  fluxes  and  the  consideration  of  additional  fluids;  (iii)  use  of  microscopy  techniques 
to  verify  contact  angles  inferred  from  experimental  data;  (iv)  meniscus  dynamic  response 
to  cyclic  heat  inputs  coupled  to  acceleration  transients;  and  (v)  meniscus  dynamic 
response  to  variations  in  the  transient  foim  of  acceleration. 

The  practical  significance  of  these  results,  as  related  to  heat  transfer  devices 
utilizing  contact  lines  and  capillarity  to  enhance  heat  and  mass  transfer,  will  be  to  provide 
information  helpful  in  defining  transient  performance  parameters  and  operational  limits. 
It  has  long  been  assumed  that  the  operating  performance  of  such  devices  subjected  to 
transverse  acceleration  components  would  not  be  altered  significantly.  However,  as 
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demonstrated  in  this  study,  the  combined  axial  and  transverse  acceleration  components 
may  result  in  an  increased  depriming  rate  leading  to  a  dry  out  condition  sooner  than 
typically  predicted  by  conventional  approaches.  Also  with  the  addition  of  heat,  the 
deprime  and  reprime  dynamics  may  be  significantly  altered  due  to  thermocapillaiy  effects 
and  result  in  either  a  degraded  or  possibly  improved  performance  such  as  with  the  arterial 
heat  pipe  discussed  by  Yeikes  et  al.^^. 
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APPENDIX  A 


MATHEMATICAL  DERIVATIONS 


Acceleration  Field 

The  acceleration  field  was  derived  using  the  geometry  shown  in  Fig.  A.  la. 

For  0  <J)  ^  , 

2 


=  zcostj) 

Tj.  =  rQ-zcos<j) 

Tj.  =  ro+  (z^-z)  cos<j) 


and 


y  =  Vq-  iz^-z)  sin<|) 


(2) 


The  position  vector  from  the  non-inertial  reference  frame 


^T^Zr^y^y 


(3) 


now  becomes 


111 


r 


(4) 


=  (ro+  {z^-z)  cos<t))^j.^+(yQ-  iz^-z)  sin<l))^y 


and  its  derivative 


di> 

dt 


-fcos* 


(-  <^1  ■  Sin4 


(5) 


with  the  angular  velocity  vector 


a  =  (oSy. 


(6) 


Referencing  an  inertial  reference  frame 


d^r 

dt2 


dfF 

dt2 


+  2(0  X 


dF 

dt 


+ 


+  (oxCcoxr^) 


(7) 


and  rearranging  for  the  non-inertial  frame  as  reference,  the  effective  forces  experienced 
by  an  observer  on  the  centrifuge  becomes 
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sff 


d}r<  d^T 

m - r-  =  Id -  +  m 


dt‘ 


dt‘ 


=  ^ 


-2(i)x-^-  -  -  tox((i)x?') 

at  dt 


(8) 


-2(0 X xr'  -  cox  ((oxr) 
dt  dt 


=  ^ 


-2(o(-^  (Zj-z)  sin(j)+-^cos<J))  +  -^(ro+ (z^-z)  cos(|)) 


dz 

dt 


d(o 

dt 


m  [(o^(ro+  (z^-z)  cos(|))]d^^ 


(9) 


If  the  capillary  tube  is  at  a  fixed  angle  of  inclination  such  that  d<t>/dt  =  0  and  assuming 
the  only  external  force  is  due  to  gravity  =  m{-g)  §y)  the  effective  force  becomes 


Fgff  =  m  [(o2(ro+  (z^-z)  cos<|))]^^^ 


+  m(.-g)  Sy  +  m 
=  mb 


d(o 


dt 


(ro+  (z^-z)  cos<J))-2(o-^cos<()  ^^0) 


Now  the  acceleration  vector  is  cast  into  a  form  referenced  to  a  non-inertial  reference 
frame  affixed  to  the  centrifuge  where 


Rotating  the  non-inertial  reference  frame  (180  -  <i))°  about  as  shown  in  Fig.  A.  lb,  the 
transformed  acceleration  vector  referenced  to  the  capillary  tube  can  be  determined  where 

F  =  Ab  (12) 
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with 


A  = 


§'  •§,  s'  S'r  'S^ 

j-  T  T  y  T  ^ 

-S  -S  •S 

^y  ^xj.  ^y  ^y  ^y 

^T’^Tr  ^T'^y  4-^r 


(13) 


A  = 


COS  (180-4)) 
cos  (90+4)) 
cos  (90) 


cos  (90-4))  cos  (90) 
cos  (180-4))  cos  (90) 
cos  (90)  cos(O) 


(14) 


A  = 


-cos4)  sin4)  0 
-sin4)  -cos4)  0 
0  0  1 


(15) 


now 


F  = 

-COs4) 

-sin4) 

sin4)  0 
-cos4)  0 

0 

0  1 

(16) 
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where 


=  (b^^{-cos^)  +byS±ni>^S'j.^ 

+  (jbj.^(-sin<j))  +jby(-cos<())  j^y  (17) 

+ 


Assuming  a  fixed  axial  location,  dz/dt  =  0,  and  eliminating  the  prime  notation  letting 

and  combining  the  dy  and  d^  components  to  form  a  magnitude  of  the  transverse 
component  djj 

b  =  [-w^  {r'o+  (z^-  z)  cos<J>)  cos<j)+  {-g)  sin4)]d^  + 

[-(i)2(ro+  (Zi-z)cos<j>)  sin<|)-  (-5-)cos4)]^  + 

(18) 

(Note,  the  assumption  of  a  fix;ed  axial  location,  dz/dt  =  0,  implies  negligible  effects 
due  to  the  Coriolis  acceleration  which  is  embodied  in  the  transverse  acceleration 
component.  The  Coriolis  acceleration  is  induced  from  the  motion  of  the  fluid  column  or 
meniscus.  Therefore  prior  knowledge  of  the  meniscus  dynamics,  specifically  tiie  velocity, 
is  required  to  evaluate  the  Coriolis  acceleration.  The  one-dimensional  equation  of  motion 
for  the  meniscus  formulated  for  this  investigation  assumes  that  the  transverse  acceleration 
component  has  a  negligible  effect  on  the  dynamic  behavior  of  the  meniscus  and  therefore 
does  not  accotmt  for  tangential  and  Coriolis  acceleration  effects.) 

If  the  angular  velocity  is  assumed  to  be  linear  as  a  function  of  time  where  the 


(ro+  {z^-  z)  cos^) 

dt 


l/c 
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transient  angular  velocity  is  of  the  form 


(0  (t)  =  Bt, 


tf 


tf 


^t^tf 


(19) 


where 


B  = 


oj 


(20) 


The  coefficient  S  can  be  determined  as  a  function  of  the  desired  peak  radial  acceleration 
2-t  t  —  izf/2 ,  cycle  time,  and  radial  position  at  the  pivot  point.  Using  the  angular 

velocity  at  t  =  can  be  found  in  terms  of  the  coefficient  B  where 


(0  = 
9pk  ~ 


(21) 


Solving  for  B  results  in 


B  = 


'■  9pk 


itf)^r 


(22) 


A  typical  angular  velocity  profile  with  a  cycle  ftequency  of  0 . 0015  is  shown 
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in  Fig.  A.2  with  a  4.2-g  peak  radial  acceleration.  The  resulting  acceleration  fields 
referenced  to  the  capillary  tube  with  the  inclination  angle  ranging  from  0"  to  60°  are 
shown  in  Fig.  A.3. 


Equation  of  Motion 

A  simple  analytical  model  was  derived  to  predict  the  motion  of  the  meniscus 
subjected  to  transient  accelerations.  The  analytical  formulation  of  the  equation  of  motion 
was  simplified  assuming  bulk  flow  only  in  the  axial  direction  of  the  capillary  tube. 
Effectively  the  liquid  column  in  the  capillary  tube  is  assumed  to  undergo  a  slug  flow. 
From  the  conservation  of  momentum, 

=  F^  +  F^  =  ,  (23) 


the  surface  and  body  forces  in  the  axial  direction  combine  resulting  in  the  axial  force 

+  (24) 

Consider  a  control  volume,  depicted  in  Fig.  A.4,  bordered  by  the  walls  of  the 
capillary  tube,  meniscus  surface,  and  a  level  of  the  surface  of  die  reservoir.  The  surface 
force,  Fg^,  is  formulated  by  accounting  for  the  net  pressure  force  on  the  meniscus  surface 

and  at  the  level  of  the  reservoir,  cumulative  effect  of  solid-liquid  intermolecular  forces 
in  the  near  contact  line  region,  and  the  wad.  shear  stress  associated  with  Poiseudle  flow 
in  the  liquid  column.  The  contact  line  force  is  defined  asstiming  a  spherical  meniscus  and 
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a  constant  equivalent  or  apparent  contact  angle  such  that 


=  o  (27ta)  cose  +  x  ^(27ta)  T]  +  (P^Tta^-PiTia^)  .  (25) 


Since  the  lower  control  volume  surface  is  at  the  level  of  the  reservoir  surface,  the 
net  pressure  force  is  zero,  =  P^.  The  shear  stress  at  the  wall  is  defined  using  the 
Darcy  friction  factor  assuming  laminar  tube  flow  where 


A 


T 


64 

Re 

4t 

w 

64  (  P  ^ave  ^ 

8“  j 

a 

4|a  dt] 
a  dt  ' 


(26) 


This  results  in  the  following  surface  force  formulation  while  using  a  slug  flow 
approximation  to  describe  the  flow  of  liquid  in  the  column,  Le.,  =  dq/dt. 

Fg  =  o  (2jca)  cosS-STcpq-^  (27) 

^  dt 
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The  total  body  force, 


(28) 


'  l''l%P(2K)rdrdz 

=  f'b  p  (na^)  dz, 


is  formulated  by  integrating  over  the  volume  of  the  liquid  column  the  axial  acceleration 
component,  Jb^,  from  eq.  (18)  as  defined  by 

=  -(0(t)2  cost})  (ro+  (Zi-z)cos<f))-g-sin<|>.  (29) 


The  transient  angular  velocity  is  of  the  form 


w  (t)  =  St, 

=  S(tf-  t)  , 


0  <;  t:^  — 
2 

2 


(30) 


where 


S  = 


^  9pJr 


4 


(  tf)  2  T, 


oj 


(31) 


Assuming  that  the  axial  fluid  velocity  is  also  defined  by  =  dt]/dt  (due  to  the  slug 

flow  assumption)  and  the  mass  in  the  capillary  tube  asi37  =  p(Txa2)Ti,the  relation  for 
the  time  rate  of  change  in  momentum  is  defined  as 
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(32) 


d{mv^) 

dt 


dr)  \2 

^  dt2  1 

^dtj 

Combining  eqs.  (27),  (28),  and  (32),  the  momentum  equation  becomes 


■  2OCOS0 

_ 

pa  _ 

pa2  '  dtj  [jo 

r  dt2  1  dtj 

Integrating  the  axial  acceleration  term, 


f^b^dz=  [-0)  ( t)  2  cos<j)(ro  +  iz^-z)  cos^)  -  srsin<J)]d2 
=  n  [w  ( t)  2  cos<|)|  ( i  )  cos4)  -  r^l  -  srsin(|) 


(34) 


rearranging,  the  momentum  equation  becomes 


2ocos6 

pa 


dt‘ 


dt 


pa‘ 


dt 


w(t)2  cos(|)j(2i--|)cos<|)  +  roJ+gsin<|) 


(35) 


with  initial  conditions 


0,  t  =  0 


(36) 


A  non-dimensional  form  of  the  momentum  equation  can  now  be  derived  using  the 
following  dimensionless  parameters 
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(37) 


c  =  0^  c  ^  1 

ho 

©  =  -^,  0^  ©  :^  1. 


The  static  wicking  height,  h^,  is  defined  using  the  static  contact  angle,  0  where 


2ocos0y 
pasrsin<|)  ’ 


(38) 


A  new  angular  velocity,  Q ,  is  formed  by  factoring  out  of  the  coefficient  B  (defined 
by  eq.  (31))  in  eq.  (30)  and  casting  it  in  terms  of  the  non-dimensional  time,  ©,  where 


^  0,  0^©:^  — 
=  f-^)2  (1-0),  -i^©<:l. 

rM  2 

\  O  / 


Q(©)  =  -A 


The  non-dimensional  form  of  the  momentum  equation  in  terms  of  a  transient 
dimensionless  meniscus  position  becomes 


2OCOS0 


pa 


dH 


tf^  d02 


.Ji-L 


8]lho  dC 

IJ 

pa^tf  d& 

Q(©) 


^pj^cos<l)||z^  -  -^jcos<|)  +  Toj  +  p-sin<|) 


(40) 


dividing  through  by  h^gpj^ 
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(41) 


2OCOS0 

K 

.dH  J  dif 

rde2  (d©)J 

1 

Q  (©)  ^cos4)| 

z.  -  ]cos<i)  +  To  1  +  -S-  sind) 

\ 

2  ;  ) 

The  momentum  equation  can  be  modified  to  include  a  constant  receding  or 
advancing  contact  angle,  0  differing  firom  that  of  the  static  contact  angle,  0  It  may 
be  noted  that  the  general  form  of  the  dynamic  contact  angle  is  known  to  be  functionally 
dependent  on  the  contact  line  receding  and  advancing  velocity,  0^  =  0j(dC/d0). 
Substituting  the  static  wicking  height,  h^,  into  the  left  hand  side  of  eq.  (41)  and  realizing 
that  the  contact  angle,  0 ,  is  equivalent  to  the  dynamic  contact  angle  where  ,  0  =  0  j, 
results  in 


2OCO30 

P  ^pk 


2ocos6, 


(  2ocos0_l 


g  COSB^ 

^  coie^ 


sin4). 


(42) 


Defining  the  following  dimensionless  niunbers 
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acceleration 
surface  tension 


Bo 


o 


viscous  1  _  PVc  Y  _ 

surface  tension  \'  o'  tf  ' 

'  inertia 
viscous 


(43) 


the  non-dimensional  fonn  of  the  momentum  equation  becomes 


g  COSOrf 

^  cose^ 


sin^ 


■  Re  Ca  ■ 

+ 

16Ca] 

k 

2Bo 

de=  1  d0/ 

Bo 

(44) 


/ 

£2  (0)  2  cos<t) 

Ui  -  )cos<|)  +  To  1  +  -^  sin4) 

V  2  ;  )  ffpk  . 

with  initial  conditions 


C  =  1, 


d9 


0,0  =  0 


(45) 


The  first  and  second  terms  of  eq.  (44)  represent  the  inertial  and  viscous  effects, 
respectively,  with  the  coefficients  formed  as  the  products  of  the  dimensionless  parameters 
of  aspect  ratio,  Reynolds,  cq)illary,  and  Bond  numbers.  The  Reynolds  and  capillary 
numbers  are  referenced  to  a  characteristic  velocity,  v^,  which  is  the  maximum  attainable 
velocity  over  one-half  of  the  acceleration  cycle  period  assuming  q  =  0  at  The 
Bond  number  is  referenced  to  the  peak  radial  acceleration  component, 

For  small  inertial  and  viscous  effects,  the  momentum  equation  can  be  reduced  to 
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the  following  quadratic  form  for  the  transient  dimensionless  meniscus  position. 


\ 

D  (0)  ^  cos<l)  i z ^cos^  ■¥  r'o)  +  ^  sin4> 

r  . 

g  COsBj  .  ) 

— - ^sin<t) 

9pk  cosB^ 

Q(e)2  _£  C0S2(() 

2  ; 

Q(0)2  i  cos2<t) 

(46) 


Discretization  was  done  using  a  central  difference  scheme  in  time  only  where 


K  - 

dd 


Ci.3-C. 

2A0 


+  O(A0)2 


(47) 


_ 

d02 


Ci+2  ^  ^i+1  C  i 


A02 


5-+O(A0)2 


(48) 


and  substituting  into  the  non-dimensional  equation  of  motion,  eq.  (43),  results  in 


COS0, 


Sp*  cos©^ 


sin^)  = 


ReCa 


2  Bo 


•i+li 


^1+2  C  i  ^  (  Ci+2  “  C  i 


A02 


16Ca] 

r  (  ^i+2 

Bo 

r 

rn  2A0  jj 

2A0 


t  J  L  \ - —  /J 

Ci+1  Q  (0)  2cos<()^^Zj^  -  — 2^ili.jcos<j)-Xoj  sin4> 


(49) 


Combining  like  terms  and  expanding  results  in 


124 


aRe 

8h^ 
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^  AQ^Bo  g  cosQg 
Ca  gp^  cose  ^ 

AQ^Bo 
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Ca 

aRe 
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h. 


aRe 


2 

i+1 


8/1, 


C 

C^i 


(50) 


with  initial  conditions 


c  =  i,  =0,  e  =  o 


(51) 


or 

Cl  =  Ci.i  =  1,  0  =  0  (52) 

Now  Cj+2  ^  determined  at  each  time  step,  A0,  by  solving  in  the  quadratic  using  the 
root  consistent  for  a  converging  solution  over  the  cycle  time,  0  ^0^1. 

Solutions  to  eq.  (50)  were  calculated  by  varying  the  Bond  and  capillary  numbers. 
Numerical  error  due  to  time  step  size  is  on  the  order  of  (A0)  ^  and  was  minimized  by 
obtaining  solutions  at  varying  time  step  size  and  selecting  a  time  step  appropriate  to  have 
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negligible  effects  on  the  resulting  solution  (A6  :s  0.0001).  These  solutions  were 
compared  to  the  inviscid  solution  and  the  percent  deviation  or  retardation  from  the 
inviscid  solution  determined.  In  these  solutions  the  static  contact  angle  of  16.7®  was 
assumed  to  be  equivalent  to  the  dynamic  contact  angle.  It  is  interesting  to  note  that  the 
requirement  that  the  static  contact  angle  be  specified  is  for  the  calculation  of  the  aspect 
ratio  and  capillary  and  Reynolds  numbers.  For  fixed  Bond  and  c^illary  numbers  and  for 
a  dynamic  contact  angle  equivalent  to  the  static  contact  angle,  solutions  are  insensitive 
to  the  changes  in  the  static  contact  angle  since  the  cosine  of  the  static  contact  angle 
cancels  out  in  the  inertial  grouping  of  the  aspect  ratio  and  Reynolds  number.  This  is  not 
the  case,  however,  when  the  dynamic  contact  angle  differs  firom  the  static  contact  angle. 
By  comparing  the  inertial  and  viscous  terms,  the  relative  magnitude  and  functional 
relationship  to  each  other  was  ascertained.  Tables  A.la-A.3a  show  the  conditions  for  the 
numerical  solutions  for  a  capillary  angle  of  inclination  of  25®,  30®,  and  50®.  Tables  A.lb- 
A.3b  show  the  relevant  statistical  information  for  logarithmic  linear  regressions  as  the 
functional  relation  between  the  inertial  and  viscous  terms  for  each  angle  of  inclination. 
Figures  A.5-A.7  graphically  present  the  functional  relationships  between  the  inertial  and 
viscous  terms  as  well  as  showing  the  magnitude  of  the  viscous  term  contributing  to  the 
dynamic  response  of  the  meniscus  resulting  in  a  retardation  of  the  meniscus  height  from 
the  inviscid  solution. 

Tables  A.4-A.7  are  the  computer  programs,  written  in  C,  to  calculate  acceleration 
fields  and  dimensionless  meniscus  position  used  for  this  investigation. 
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a) 


Fig.  A.l  Capillary  tube  orientation:  a)  pivot  point  location,  directional  coordinates  and 
b)  coordinate  systems  as  referenced  to  the  centrifuge  accelerometer  and  capillary  tube. 
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Fig.  A.2  Plot  of  the  centrifuge  angular  velocity  test  cycle  {f  =  0.0015  Hz). 
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Axial  and  Transverse  Acceleration  Components 
20°  inclination,  f  =  0.001  5  Hz,  gpk  =  4.2  g 


Axial  and  Transverse  Acceleration  Components 

40°  inclination,  f  =  0.001  5  Hz,  gpk  =  4.2  g 


bR  (magnitude) 


O 


Fig.  A.3(contd)  c)  =  40°, 
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Axial  and  Transverse  Acceleration  Components 
60°  inclination,  f  =  0.001  5  Hz,  gpk  =  4.2  g 


bR  (magnitude) 


Fig,  A.3(contd)  d)  <j>  =  60°, 


132 


Table  A.la 

Dimensionless  parameters  and  %  deviation 
from  the  inviscid  solution  (<|)  =  25  ®). 
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Table  A.lb 

Statistical  infoimation  showing  the  functional  relationship 
between  inertial  and  viscous  terms  (<})  =  25  °). 


Bo 

4> 

ind.  var. 

dep.  var. 

R-sqd. 

slope 

intercept 

la 

log 

a  i?eCa'\ 

log(^) 

1.0 

0.50 

-0.26007 

0.4 

(K  BO  j 

\Boj 

1.0 

-0.63634 

1.0 

1.0 

-1.1338 

1.5 

1.0 

-1.3539 

aU 

25" 

intercepts 

log (So) 

1.0 

-1.25 

-1.1338 
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viscous  term,  (Ca/Bo) 


— 6H-  Bo  =  0.2 
— e —  Bo  =  0.4 
—>1^  Bo  =  1  .0 
— A —  Bo  =  1 .5 

log(Ca/Bo)  =  -(1 .1  338  +  1 .25*log(Bo))  +  0.50*log((Re*Ca*a)/(Bo*ho)) 

]  Q-2  t — I  "r rnTTTi 1  i  i  i  vttti — i  i  i  i  iftfi 1  i  i  mii| — ttttttttt 1  i  i  i  nni  i  i  1 1  imi  r-rnTn^ 

:  viscous  terms  ^ 

-  contributing  to  the  dynamic 

1 0~3  ^  response  of  meniscus’  1 


1  0-4 


1  0-5 


1  0-6 


1  0-7 


*  >0.7%  deviation  from  inviscid  solution 
at  end  of  cycle  for  25°  inclinotion 


0-11  10-10  10-9  10-8  10-7  10-6  10-5  10-4  10-5 

inertial  term,  (Re*Ca*a)/(Bo*ho) 


Fig.  A.5  Functional  relationship  between  viscous  and  inertial  terms  (<|)  =  25  °). 


Table  A.2a 

Dimensionless  parameters  and  %  deviation 
from  the  inviscid  solution  (<|)  =  30  °). 


<J> 

a  ReCa 

K  Bo 

Ca 

Bo 

%  deviation 

Em 

o 

0 

9.797e-5 

5.0e-3 

41.4 

l.Oe-4 

9.797e-7 

5.0e-4 

8.9 

mm 

9.797e-9 

5.0e-5 

0.24 

l.Oe-6 

9.797e-ll 

5.0e-6 

0 

0.4 

l.Oe-3 

1.386e-4 

2.5e-3 

30.5 

l.Oe-4 

1.386e-6 

2.5e-4 

3.8 

l.Oe-5 

1.386e-8 

2.5e-5 

0.06 

l.Oe-6 

1.386e-10 

2.5e-6 

0 

1.0 

l.Oe-3 

2.191e-4 

l.Oe-3 

16.9 

l.Oe-4 

2.191e-6 

l.Oe-4 

0.86 

l.Oe-5 

2.191e-8 

l.Oe-5 

0.01 

l.Oe-6 

2.191e-10 

l.Oe-6 

0 

1.5 

l.Oe-3 

2.683e-4 

6.7e-4 

11.7 

l.Oe-4 

2.683e-6 

6.7e-5 

0.37 

l.Oe-5 

2.683e-8 

6.7e-6 

0 

l.Oe-6 

2.683e-10 

6.7e-7 

0 
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Table  A.2b 

Statistical  infonnation  showing  the  functional  relationship 
between  inertial  and  viscous  terms  (<})  =  30  °). 


Bo 

4) 

ind.  var. 

dep.  var. 

R-sqd. 

slope 

intercept 

0.2 

1  nrr 

log(-^) 

1.0 

0.50 

-0.29658 

0.4 

\Bol 

1.0 

-0.67294 

1.0 

1.0 

-1.1703 

1.5 

1.0 

-1.3904 

all 

30° 

intercepts 

log  (Bo) 

1.0 

-1.25 

-1.1703 
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viscous  term,  (Ca/Bo) 


Table  A.3a 

Dimensionless  parameters  and  %  deviation 
from  the  inviscid  solution  (<|)  =  50  ®). 


Bo 

Ca 

B 

a  ReCa 

K  Bo 

Ca 

Bo 

%  deviation 

m 

1.501e-4 

5.0e-3 

29.0 

DSH 

1.501e-6 

5.0e-4 

3.3 

UsiSH 

1.501e-8 

5.0e-5 

0.05 

l.Oe-6 

1.501e-10 

5.0e-6 

0 

n 

l.Oe-3 

2.123e-4 

2.5e-3 

18.7 

LOe-4 

2.123e-6 

2.5e-4 

1.1 

2.123e-8 

2.5e-5 

0.01 

l.Oe-6 

2.123e-10 

0 

1.0 

l.Oe-3 

3.356e-4 

l.Oe-3 

7.8 

l.Oe-4 

3.356e-6 

l.Oe-4 

0.18 

l.Oe-5 

3.356e-8 

0 

l.Oe-6 

3.356e-10 

l.Oe-6 

0 

1.5 

l.Oe-3 

4.111e-4 

6.7e-4 

4.3 

4.111e-6 

6.7e-5 

0.06 

l.Oe-5 

4.111e-8 

6.7e-6 

0 

l.Oe-6 

4.111e-10 

6.7e-7 

0 

- — - 
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Table  A.3b 

Statistical  information  showing  the  functional  relationship 
between  inertial  and  viscous  terms  (<j)  =  50  °). 


Bo 

4> 

ind.  var. 

dep.  var. 

R-sqd. 

slope 

intercept 

0.2 

50° 

loaf  ^ 

log(-2a) 

1.0 

0.50 

-0.38922 

0.4 

so  ) 

\Bol 

1.0 

-0.76553 

1.0 

1.0 

-1.2629 

1.5 

1.0 

aU 

50° 

intercepts 

log  (Bo) 

1.0 

-1.25 

-1.2629 
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iscous  term,  (Co/Bo) 


> 


Bo  =  0.4 
Bo  =  1 .0 
Bo  =  1 .5 


log(Ca/Bo)  =  -(1 .2629  +  1 .25*log(Bo))  +  0.50»log((Re*Ca*a)/(Bo*ho)) 


10-11  10-10  10-9  1 0-8  10-7  10-6  10-5  1 0-'i  10-5 


inertial  term,  (Re*Ca*a)/(Bo*ho) 


Fig.  A.7  Functional  relationship  between  viscous  and  inertial  terms  (<|)  =  50  ®). 


Table  A.4 

Computer  program  to  calculate  the  acceleration  field  in  space  and  time. 


#include  <stdio.h> 
#include  <math.h> 


main() 

{ 


float  angularfireq; 

/♦angular  velocity  frequency,  Hz*/ 

float  tf; 

/♦time  period  for  angular  velocity  cycle,  sec*/ 

float  accelpk; 

/♦peak  acceleration  at  tf/2,  "g"*/ 

float  ro; 

/♦radial  location  of  cap  tube  pivot  point,  m*/ 

float  rt; 

/♦radial  table  coordinate*/ 

float  angularvell; 

/♦angular  velocity  for  t  =  O-tf/2*/ 

float  angularvel2; 

/♦angular  velocity  for  t  =  tf/2-tf*/ 

float  t; 

/♦time,  sec*/ 

float  B; 

/♦positive  slope  of  angular  velocity*/ 

float  zl; 

/*zl  =  z/L  =  0.5*/ 

float  z; 

/♦z*L  =  axial  cap  tube  coordinate,  m*/ 

float  angle; 

/♦inclination  angle  of  cap  tube,  radians*/ 

float  L; 

/♦cap  tube  length,  m*/ 

float  bz; 

/♦axial  acceleration  component,  "g"*/ 

float  bR; 

/♦transverse  acceleration  component,  "g"*/ 

float  k,i; 

char  filenm[20]; 

FILE 

piintf  ("\nThis  program  wiU  calculate  the  axial  and  transverseNn"); 
printf  ("acceleration  components  as  a  function  of  q)ace  and  timeNn"); 
printf  (’Venter  inclination  angle,  deg.'^"); 
scanf  ("%f’,&angle); 
angle  =  angle*3. 141592/180; 

printf  ("\nenter  cycle  frequency  for  angular  velocity,  HzV); 
scanf  ("%f’,&angularfreq); 

printf  ("\nenter  peak  acceleration  (half  period),  g’sNn"); 

scanf  ("%f',&accelpk); 

printf  ("\nenter  filename  for  your  dataNn"); 

scanf  ("%s",filenm); 
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^  =  fopen  (filemn,”w"); 
ro  =  0.8073; 

L  =  10/100; 
tf  =  1/angularfreq; 

B  =  sqrt(accelpk/(ro*pow(tf/2,2))); 

zl  =  0.5; 
t  =  0; 

for  (i=0;  i<=tf;  i++) 

{ 

z  =  0; 


for  (k=0;  k<=10;  k++) 

{ 

rt  =  ro  +  (zl-z)*L*cos(angle); 

/*(t  =  O-tf/2)*/ 
if  (t  <=  tf/2) 

{ 

angularvell  =  B*t; 

bz  =  -  pow(anguiarvell,2)*rt*cos(angle)  +  (-9.81)*sin(angle)/9.81; 

Br  =  pow(pow(-  pow(anguIarvell,2)*rt*sin(angle)  -  (-9.81)*cos(angle)/9  81  2) 
+  pow(rt*B,2),0.5); 

} 


/*(t  =  tf/2-tf)*/ 
if  (t  >  tf/2) 

{ 

angularvel2  =  B*(tf-t); 

bz  =  -  pow(angularvel2,2)*it*cos(angle)  +  (-9.81)*sin(angle)/9.81; 

Br  =  pow(pow(-  pow(angularvel2,2)*rt*sin(angle)  -  (-9.81)*cos(ailgle)/9  81  2) 
+  pow(rt*(-B),2),0.5); 

} 

^rintf  (^,"%4.3eNt%4.3eNt%4.3e\t%4.3e\n",t,z,bz3r); 

z  =  z  +  0.1; 

} 

t  =  t  +  1; 

} 

fclose  (^); 

} 
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Table  A.5 

Computer  program  to  calculate  the  acceleration  at  the  capillary 
tube  pivot  point  in  time. 


#include  <stdio.h> 
#include  <math.h> 


main() 

{ 


float  angularfreq; 
float  tf; 
float  accelpk; 
float  ro; 
float  rt; 

float  angularvell; 

float  angularvel2; 

float  t; 

float  B; 

float  zl; 

float  z; 

float  angle; 

float  L; 

float  bz; 

float  Br; 

float  dt; 

float  k,i; 


/*angular  velocity  frequency,  Hz*/ 

/*time  period  for  angular  velocity  cycle,  sec*/ 
/*peak  acceleration  at  tf/2,  "g"*/ 

/*radial  location  of  cap  tube  pivot  point,  m*/ 
/*radial  table  coordinate*/ 

/*angular  velocity  for  t  =  O-tf/2*/ 

/*angular  velocity  for  t  =  tf/2-tf*/ 

/*time,  sec*/ 

/*positive  slope  of  angular  velocity*/ 

/*zl  =  z/L  =  0.5*/ 

/*z*L  =  axial  c^  tube  coordinate,  m*/ 
/*inclination  an^e  of  cap  tube,  radians*/ 

/*cap  tube  length,  m*/ 

/*axial  acceleration  component,  "g"*/ 
/♦transverse  acceleration  component,  "g"*/ 
/♦time  step*/ 


char  filenm[20]; 
FILE  *fp; 


printf  (’^Nthis  program  will  calculate  the  axial  and  transverseNn"); 

printf  ("acceleration  components  as  a  function  of  space  and  time\n"); 

printf  ("Nnenter  inclination  angle,  deg.'tfi"); 

scanf  ("%f ',&angle); 

angle  =  angle*3. 141592/180; 

printf  (’Nnenter  cycle  frequency  for  angular  velocity,  Hz^"); 
scanf  ("%f',&angularfreq); 

printf  ('Venter  peak  acceleration  (half  period),  g’sNn"); 
scanf  ("%f ',&accelpk); 
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printf  ("\nenter  time  stepNn"); 
scanf  ("%f ’,&dt); 

printf  ("\nenter  filename  for  your  dataNn"); 
scanf  ("%s",filenm); 

^  =  fopen  (filenm,"w"); 
ro  =  0.8073; 

L  =  10/100; 
tf  =  1/angulaifreq; 

B  =  sqrt(accelpk/(ro*pow(tf/2,2))); 
zl  =  0.5; 


t  =  0; 
do 
{ 

z  =  zl; 


rt  =  ro  +  (zl-z)*L*cos(angle); 

/*(t  =  O-tf/2)*/ 
if  (t  <=  tf/2) 

{ 

angularvell  =  B*t; 

bz  =  -  pow(angularvell,2)*rt*cos(angle)  +  (-9.81)*sin(angle)/9.81; 

Br  =  pow(pow(-  pow(angiilarvell,2)*rt*sin(angle)  -  (-9.81)*cos(angle)/9.81,2) 
+  pow(rt*B,2),0.5); 

} 


/*(t  =  tf/2-tf)*/ 
if  (t  >  tfi2) 

{ 

angularvel2  =  B*(tf-t); 

bz  =  -  pow(angularvel2,2)*rt*cos(angle)  +  (-9.81)*sin(angle)/9.81; 

Br  =  pow(pow(-  pow(angularvel2,2)*rt*sin(angle)  -  (-9.81)*cos(angle)/9.81,2) 
+  pow(rt*(-B),2),0.5); 

} 

fprintf  (fp,"%4.3e\t%4.3eWo4.3eNt%4.3eV’,t,z,bz3r); 
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t  =  t  +  dt; 

} 

while  (t<=tf); 
fclose  (^); 

} 


Table  A.6 

Computer  program  to  calculate  the  transient 
dimensionless  meniscus  position  for  ethyl  alcohol. 


#include  <stdio.h> 
#include  <math.h> 


main() 

{ 


double  angularfreq;  /*angular  velocity  frequency,  Hz*/ 


double  tf; 
double  accelpk; 
double  ro; 
double  rt,rtv; 
double  t; 
double  B; 
double  zl; 
double  z; 
double  angle; 
double  dangle; 
double  thetas; 
double  theta; 

double  thetaa; 

double  L; 
double  a; 
double  d; 
double  gamma; 
double  rho; 
double  vise; 
double  ho; 
double  Vc; 
double  Bo; 
double  Ca; 
double  Re; 
double  H0,Hl,H0v,Hlv; 
double  Hinert; 


/♦time  period  for  angular  velocity  cycle,  sec*/ 
/♦peak  acceleration  at  tf/2,  "g"*/ 
/♦radial  location  of  cap  tube  pivot  point,  m*/ 
/♦radial  table  coordinate*/ 

/♦time,  sec*/ 

/♦positive  slope  of  angular  velocity*/ 

/♦zl  =  z/L  =  0.5*/ 

/*z*L  =  axial  cap  tube  coordinate,  m*/ 
/♦inclination  angle  of  cap  tube,  radians*/ 
/♦inclination  angle  increment*/ 

/♦static  contact  angle,  radians*/ 
/♦apparent  recession  contact  angle 
(assumed  constant),  radians*/ 

/♦apparent  advancing  contact  angle 
(assumed  constant),  radiam*/ 

/♦cap  tube  length,  m*/ 

/♦cap  tube  radius,  m*/ 

/♦cap  tube  diameter,  m*/ 

/♦surface  tension,  N/m*/ 

/♦fluid  density,  kg/cu.m*/ 

/♦fluid  viscosity,  kg/m-s*/ 

/♦static  wicking  height,  m*/ 

/♦characteristic  velocity,  m/s*/ 

/♦Bond  number*/ 

/♦Capillary  number*/ 

/♦Reynolds  number*/ 

/♦dimensionaless  height,  h/ho*/ 
/♦dimensionaless  height  inviscid 
with  inertia  term,  h/ho*/ 
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double  Hvisc; 


/♦dimensionaless  height  with 
viscous  and  inertia  terms,  h/ho*/ 
double  Happrs;  /*dimensionaless  height  approximated 

for  small  inertia  and  viscous  terms,  h/ho*/ 
double  acoef,bcoef,ccoef;/*coefficients*/ 
double  al,bl,cl;  /*coef¥icients*/ 
double  alv,blv,clv;  /*coefficients*/ 
double  T;  /*dimensionaless  time,  i/tf*/ 

double  time;  /*dimensionaless  time,  t/tf  (muli^le  cycles)*/ 

double  Dt;  /*incremental  dimensionaless  time*/ 

double  omegal;  /*dimensionaless  angular 

velocity,  T  =  0-0.5*/ 

double  omega2;  /*dimensionaless  angular 

velocity,  T  =  0.5-1. 0*/ 

float  cyclenum;  /*number  of  cycles*/ 

float  k,i,j,n; 

char  £ilenm[20]; 

FILE  *fp; 

printf  ('\Nthis  program  will  calculate  the  Nn"); 
printf  ("meniscus  height  as  a  function  of  timeSn"); 
printf  ("\nBond  number  (eg.  0.1-2.0)  =  "); 
scanf  ("%lf',&Bo); 

printf  ("NnCapfllary  number  (eg.  le-4  -  le-8)  =  ”)• 
scanf  ("%lf' ,&Ca); 

printf  ("Snenter  peak  acceleration  (half  period),  g’sNn"); 
scanf  ("%lf',&accelpk); 

Bo  =  Bo/accelpk; 

printf  ('^ente^  number  of  cyclesVi"); 
scanf  ("%lf' ,&cyclenum); 

printf  ("Nnenter  dimensionaless  incremental  time  step,  Dt  <  0.00  INn")- 
scanf  ("%lf ',&Dt); 

printf  ("Nnenter  inclination  angle,  deg.Nn"); 
scanf  ("%lf',&angle); 
angle  =  angle*3.141592/180; 
printf  (’\ienter  static  contact  angle,  deg.'^"); 
scanf  ("%lf',&thetas); 

printf  ("Center  recession  contact  angle,  deg.Nn"); 
scanf  ("%lf' ,&theta); 

printf  ("Nnenter  advancing  contact  angle,  deg.Vi"); 
scanf  ("%lf ',&thetaa); 
thetas  =  thetas*3. 141592/1 80; 
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theta  =  theta*3.141592/180; 
thetaa  =  thetaa*3. 141592/180; 
printf  ("\nenter  filename  for  your  dataNn"); 
scanf  ("%s",filenm); 

=  fopen  (filenm,"w"); 
ro  =  0.8073; 

L  =  10/100; 
zl  =  0.5; 

rho  =  781.0;  /*Ethyl  alcohol  at  30C*/ 
gamma  =  .02189;  /*Ethyl  alcohol  at  30C*/ 
vise  =  0.001;  /*Ethyl  alcohol  at  30C*/ 
a  =  sqrt(Bo*gamma/(4*rho*9.81)); 
d  =  2*a; 

Vc  =  gamma*Ca/visc; 

ho  =  2*gamma*cos(thetas)/(rho*a*9.81*sin(angle)); 
tf  =  2*ho/Vc; 

Re  =  rho*Vc*d/visc; 


B  =  sqrt(accelpk/(ro*pow(tf/2,2))); 


i^rintf 

(^,"NntimeNtHvisc\tHinert\tHapprx>tBoNtCa''lReNta/hdStangleNtinertiaNtviscousNn"); 


/*initial  conditions*/ 
HOv  =  1.0; 

Hlv  =  HOv; 

HO  =  1.0; 

HI  =  HO; 

Happrx  =  HO; 

T  =  0; 
time  =  0; 
k  =  0; 

n  =  cyclenum; 


printf 
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(^;'%43eNt%43eWo43e\t%4.3eNt%4.3e\t%4.3e'«%4.3e\t%4.3eNt%4.3e'«%4.3eNt%4.3e^n", 

time  ,HOv, HO, Happrx,Bo*accelpk,Ca, Re, a/ho, angle/3. 141592*180,Re*Ca*a/Bo/accelpk/ 
ho,Ca/Bo/accelpk); 

do 

{ 


for  (i=0;  i<=l/Dt;  i-H-) 

{ 


T  =  T  +  Dt; 
k  =  k  +  Eh; 
time  =  time  +  Dt; 

It  =  ro  +  (zl*L  -  ho*Hl/2)*cos(angle); 
rtv  =  ro  +  (zl*L  -  ho*Hlv^)*cos(angle); 

/*(T  =  0-0.5)*/ 
if  (T  <=  0.5) 

{ 

omegal  =  B*tf*T; 

/*viscous  foimuladon  with  inertia  term*/ 
alv  =  a*Re/(8*ho); 

blv  =  (a*Re*Hlv/(2*ho))  +  (8*Dt*Hlv)  -  (a*Re*H0v/(4*ho)); 

civ  =  -(Dt*Dt*Bo/Ca  *  (sin(angle)*cos(theta)/cos(thetas) 

-  Hlv*(pow(omegal,2)*cos(angle)*rtv  +  sin(angle))) 

-  (a*Re/(2*ho)  -  8*Dt)  *  HOv  *  Hlv 

-  (a*Re/(8*ho))  *  HOv  *  HOv 
+  (a*Re/(ho))  *  Hlv  *  Hlv); 


Hvisc  =  (-blv  +  sqrt(pow(blv,2)-4*alv*clv))/(2*alv); 
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/*inviscid  formulation  with  inertia  term 


bl  =  4*H1  -  2*H0; 

cl  =  -(8*Dt*Dt*ho*Bo/(a*Re*Ca)*(sin(angle)*cos(theta)/cos(thetas) 
-  Hl*(pow(omegal,2)*rt*cos(angle)  +  sin(angle))) 

+(8*H1*H1  -  4*H1*H0  -  H0*H0)); 


Hinert  =  (-bl  +  sqrt(pow(bl,2)-4*cl))/2;  */ 


/♦approximation  without  viscous  and  inertia  terms*/ 

bcoef  =  -  (2*zl*L/ho  +  2*ro/(ho*cos(angle)) 

+  2*sin(angle)/(pow(omegal  ,2)*pow(cos(angle),2)*ho)); 
ccoef  =  2*sin(angle)*cos(theta)/cos(thetas)/(pow(omegal  ,2)*pow(cos(angle),2)*ho); 

Happrx  =  (-bcoef  -  sqrt(pow(bcoef,2)-4*ccoef))/2; 

} 


/*(T  =  0.5-1. 0)*/ 
if  (T  >  0.5) 

{ 

omega2  =  B*tf*(l-T); 

/♦viscous  formulation  with  inertia  term*/ 
alv  =  a*Re/(8*ho); 

blv  =  (a*Re*Hlv/(2*ho))  +  (8*Dt*Hlv)  -  (a*Re*H0v/(4*ho)); 

civ  =  -(Dt*Dt*Bo/Ca  *  (sin(angle)*cos(thetaa)/cos(thetas) 

-  Hlv*(pow(omega2,2)*cos(angle)*rtv  +  sin(angle))) 

-  (a*Re/(2*ho)  -  8*Dt)  *  HOv  *  Hlv 

-  (a*Re/(8*ho))  *  HOv  *  HOv 
+  (a*Re/(ho))  *  Hlv  *  Hlv); 
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Hvisc  =  (-blv  +  sqrt(pow(blv,2)-4*alv*clv))/(2*alv); 


/*inviscid  formulation  with  inertia  term 


bl  =  4*H1  -  2*H0; 


cl  =  -(8*Dt*Dt*ho*Bo/(a*Re*Ca)*(sin(angle)*cos(thetaa)/cos(thetas) 
-  Hl*(pow(omega2,2)*rt*cos(angle)  +  sin(angle))) 

+(8*H1*H1  -  4*H1*H0  -  H0*H0)); 


Hinert  =  (-bl  +  sqrt(pow(bl,2)-4*cl))/2;  */ 


/♦approximation  without  viscous  and  inertia  terms*/ 

bcoef  =  -  (2*zl*L/ho  +  2*ro/(ho*cos(angle)) 

+  2*sin(angle)/(pow(omega2,2)*pow(cos(angle),2)*ho)); 

ccoef  = 

2*sin(angle)*cos(thetaa)/cos(thetas)/(pow(omega2,2)*pow(cos(angle),2)*ho); 

Happrx  =  (-bcoef  -  sqrt(pow(bcoef,2)-4*ccoef))/2; 

} 

if  (k  >=  0.001) 

{ 

i^rintf 

(^,"%4.3eNt%4.3e\t%4.3e^t%4.3eNt%4.3eNt%4.3eNt%4.3eW<4.3eNt%4.3eNt%4.3e\t%4.3eNn", 

time, Hvisc, Hinert, Happrx3o*accelpk,Ca, Re, a/ho, angle/3.141592*180,Re*Ca*a/ho/Bo/a 
ccelpk,Ca/Bo/accelpk); 
k  =  0; 

} 

HOv  =  Hlv; 

Hlv  =  Hvisc; 

HO  =  HI; 

HI  =  Hinert; 
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} 

while  (n>0); 

fclose  (i^); 

} 


Table  A.7 

Computer  program  to  calculate  the  transient 
dimensionless  meniscus  position  for  water. 


#include  <stdio.h> 
#include  <math.h> 


main() 

{ 


double  angularfreq;  /*angular  velocity  frequency,  Hz*/ 

double  tf;  /*time  period  for  angular  velocity  cycle,  sec*/ 

double  accelpk;  /*peak  acceleration  at  tf/2,  "g"*/ 

double  ro;  /*radial  location  of  c^  tube  pivot  point,  m*/ 

double  rt^tv;  /*radial  table  coordinate*/ 

double  t;  /*time,  sec*/ 

double  B;  /*positive  slope  of  angular  velocity*/ 

double  zl;  /*zl  =  z/L  =  0.5*/ 

double  z;  /*z*L  =  axial  c^  tube  coordinate,  m*/ 

double  angle;  /*mclination  angle  of  cap  tube,  radians*/ 

double  dangle;  /*inclination  angle  increment*/ 

double  thetas;  /*static  contact  angle,  radians*/ 

double  theta;  /*jq)parent  recession  contact  angle 

(assumed  constant),  radians*/ 

double  thetaa;  /*apparent  advancing  contact  angle 

(assumed  constant),  radians*/ 
double  L;  /*cap  tube  length,  m*/ 

double  a;  /*cap  tube  radius,  m*/ 

double  d;  /*cap  tube  diameter,  m*/ 

double  gamma;  /*surface  tension,  N/m*/ 

double  rho;  /*fluid  density,  kg/cu.m*/ 

double  vise;  /*fluid  viscosity,  kg/m-s*/ 

double  ho;  /*static  wicking  height,  m*/ 

double  Vc;  /*characteristLc  velocity,  m/s*/ 

double  Bo;  /*Bond  number*/ 

double  Ca;  /*Capillary  number*/ 

double  Re;  /*Reynolds  number*/ 

double  H0,Hl,H0v,Hlv;  /*dimensionaless  height,  h/ho*/ 

double  Hinert;  /*dimensionaless  height  inviscid 

with  inertia  term,  h/ho*/ 
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double  Hvisc; 


/♦dimensionaless  height  with 
viscous  and  inertia  terms,  h/ho*/ 
double  Happrx;  /*dimensionaless  height  approximated 

for  small  inertia  and  viscous  terms,  h/ho*/ 
double  acoef,bcoef,ccoef;/*coefficients*/ 
double  al,bl,cl;  /*coefficients*/ 
double  alv,blv,clv;  /*coefficients*/ 
double  T;  /*dimensionaless  time,  t/tf*/ 

double  time;  /*dimensionaless  time,  t/tf  (mulitple  cycles)*/ 

double  Dt;  /*incremental  dimensionaless  time*/ 

double  omegal;  /*dimensionaless  angular 

velocity,  T  =  0-0.5*/ 

double  omega2;  /*dimensionaless  angular 

velocity,  T  =  0.5-1. 0*/ 

float  cyclenum;  /*number  of  cycles*/ 

float  k,i,j,n; 

char  filenm[20]; 

FILE  *^; 

printf  ("\Nthis  program  will  calculate  the  Nn"); 
printf  ("meniscus  height  as  a  function  of  timeNn"); 
printf  ('NNbond  number  (eg.  0.1-2.0)  =  "); 
scanf  ("%lf’,&Bo); 

printf  ("\Ncapillary  number  (eg.  le-4  -  le-8)  =  "); 
scanf  ("%lf’ ,&Ca); 

printf  ('Nnenter  peak  acceleration  (half  period),  g’sNn"); 
scanf  ("%lf ’,&accelpk); 

Bo  =  Bo/accelpk; 

printf  ('Vnter  number  of  cyclesVi"); 
scanf  ("%lf’,&cyclenum); 

printf  ("Nnenter  dimensionaless  incremental  time  step,  Dt  <  O.OOlNn"); 
scanf  (”%lf ',&Dt); 

printf  ("Nnenter  inclination  angle,  deg.'^"); 
scanf  ("%lf' ,&angle); 
angle  =  angle*3.141592/180; 
printf  ("Nnenter  static  contact  angle,  deg.Vi"); 
scanf  ("%lf',&thetas); 

printf  ("Nnenter  recession  contact  angle,  deg.Nn"); 
scanf  ("%lf ',&theta); 

printf  ("Nnenter  advancing  contact  angle,  deg.Nn"); 
scanf  ("%lf' ,&thetaa); 
thetas  =  thetas*3.141592/180; 


theta  =  theta*3.141592/180; 

thetaa  =  thetaa*3. 141592/180; 

printf  ("\nenter  filename  for  your  data\n"); 

scanf  ("%s",filenm); 

fp  =  fopen  (filenm,"w"); 

ro  =  0.8073; 

L  =  10/100; 
zl  =  0.5; 

rho  =  997.0;  /*water  at  300K*/ 
gamma  =  .0717;  /*water  at  300K*/ 
vise  =  855e-6;  /*water  at  300K*/ 
a  =  sqrt(Bo*gamma/(4*rho*9.81)); 
d  =  2*a; 

Vc  =  gamma*Ca/visc; 

ho  =  2*gamma*cos(thetas)/(rho*a*9.81*sin(angle)); 
tf  =  2*ho/Vc; 

Re  =  iho*Vc*d/visc; 


B  =  sqrt(accelpk/(ro*pow(tf/2,2))); 


fprintf 

(^/\itime\Thvisc\Thinert\Thapprx\Tbo\Tca\TreNta/hdNtangleNtinertiaNtviscous^''); 


/♦initial  conditions*/ 
HOv  =  1.0; 

Hlv  =  HOv; 

HO  =  1.0; 

HI  =  HO; 

Happrx  =  HO; 

T  =  0; 
time  =  0; 
k  =  0; 

n  =  cyclenum; 


l^rintf 
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(^,"%4.3e'«%4.3eWo4.3eNt%4.3e^%4.3e\t%4.3e'\t%4.3eNt%4.3eWo4.3eNt%4.3e\t%4.3eNn", 

time,H0v,H0,Happrx,Bo*accelpk,Ca,Re,a/ho,angle/3.141592*180,Re*Ca*a/Bo/accelpk/ 

ho,Ca/Bo/accelpk); 

do 

{ 


for  (i=0;  i<=l/Dt;  i-H-) 

{ 


T  =  T  +  Dt; 
k  =  k  +  Dt; 
time  =  time  +  Dt; 

rt  =  ro  +  (zl*L  -  ho*Hl/2)*cos(angle); 
rtv  =  ro  +  (zl*L  -  ho*Hlv^)*cos(angle); 

/*(T  =  0-0.5)*/ 
if  (T  <=  0.5) 

{ 

omegal  =  B*tf*T; 

/♦viscous  formulation  with  inertia  term*/ 
alv  =  a*Re/(8*ho); 

blv  =  (a*Re*Hlv/(2*ho))  +  (8*Dt*Hlv)  -  (a*Re*H0v/(4*ho)); 

civ  =  -(Dt*Dt*Bo/Ca  *  (sin(angle)*cos(theta)/cos(thetas) 

-  Hlv*(pow(omegal,2)*cos(angle)*rtv  +  sin(angle))) 

-  (a*Re/(2*ho)  -  8*Dt)  *  HOv  *  Hlv 

-  (a*Re/(8*ho))  *  HOv  *  HOv 
+  (a*Re/(ho))  *  Hlv  *  Hlv); 


Hvisc  =  (-blv  +  sqrt(pow(blv,2)-4*alv*clv))/(2*alv); 
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/*inviscid  fonnulation  with  inertia  term 


bl  =  4*H1  -  2*H0; 


cl  =  -(8*Dt*Dt*ho*Bo/(a*Re*Ca)*(sin(angle)*cos(theta)/cos(thetas) 
-  Hl*(pow(omegal,2)*rt*cos(angle)  +  sin(angle))) 

+(8*H1*H1  -  4*H1*H0  -  H0*H0)); 


Hinert  =  (-bl  +  sqrt(pow(bl,2)-4*cl))/2;  */ 


/♦^proximation  without  viscous  and  inertia  terms*/ 

bcoef  =  -  (2*zl*L/ho  +  2*ro/(ho*cos(angle)) 

+  2*sin(angle)/(pow(omegal  ,2)*pow(cos(angle),2)*ho)); 
ccoef  =  2*sin(angle)*cos(theta)/cos(thetas)/(pow(omegal  ,2)*pow(cos(angle),2)*ho); 

Happrx  =  (-bcoef  -  sqrt(pow(bcoef,2)-4*ccoef))^; 

) 


/*(T  =  0.5-1. 0)*/ 
if  (T  >  0.5) 

{ 

omega2  =  B*tf*(l-T); 

/♦viscous  formulation  with  inertia  term*/ 
alv  =  a*Re/(8*ho); 

blv  =  (a*Re*Hlv/(2*ho))  +  (8*Dt*Hlv)  -  (a*Re*H0v/(4*ho)); 

civ  =  -(Dt*Dt*Bo/Ca  *  (sin(angle)*cos(thetaa)/cos(thetas) 

-  Hlv*(pow(omega2,2)*cos(angle)*rtv  +  sin(angle))) 

-  (a*Re/(2*ho)  -  8*Dt)  *  HOv  *  Hlv 

-  (a*Re/(8*ho))  *  HOv  *  HOv 
+  (a*Re/(ho))  *  Hlv  *  Hlv); 
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Hvisc  =  (-blv  +  sqrt(pow(blv,2)-4*alv*clv))/(2*alv); 


/♦inviscid  formulation  with  inertia  term 
bl  =  4*H1  -  2*H0; 

cl  =  -(8*Dt*Dt*ho*Bo/(a*Re*Ca)*(sin(angle)*cos(thetaa)/cos(thetas) 
-  Hl*(pow(omega2,2)*rt*cos(angle)  +  sin(angle))) 

+(8*H1*H1  -  4*H1*H0  -  H0*H0)); 


Hinert  =  (-bl  +  sqrt(pow(bl,2)-4*cl))/2;  */ 


/♦approximation  without  viscous  and  inertia  terms*/ 

bcoef  =  -  (2*zl*L/ho  +  2*ro/(ho*cos(angle)) 

+  2*sin(angle)/(pow(omega2,2)*pow(cos(angle),2)*ho)); 

ccoef  = 

2*sin(angle)*cos(thetaa)/cos(thetas)/(pow(omega2,2)*pow(cos(angle),2)*ho); 

Happrx  =  (-bcoef  -  sqrt(pow(bcoef,2)-4*ccoef))/2; 

} 

if  (k  >=  0.001) 

{ 

fprintf 

(Q),"%4.3eNt%4.3eNt%4.3e\t%4.3eNt%4.3e'^%4.3eNt%4.3eNt%4.3eNt%4.3eNt%4.3eNt%4.3eSn", 

time,Hvisc,Hinert,Happtx3o*accelpk,Ca,Re,a/ho, angle/3 . 141592* 1 803e*Ca*a/ho/Bo/a 
ccelpk,Ca/Bo/accelpk); 
k  =  0; 

} 

HOv  =  Hlv; 

Hlv  =  Hvisc; 

HO  =  HI; 

HI  =  Hinert; 


160 


} 

while  (n>0); 


fclose  (fp); 

} 


APPENDIX  B 


EXPERIMENTAL  SETUP  AND  INSTRUMENTATION 

Centrifuge  Test  Bed 

Transient  acceleration  fields  were  generated  using  a  2.4-m  diameter  centrifuge 
table  (as  discussed  in  chapters  I-III)  rotating  with  a  time  variant  angular  velocity.  Single 
and  multiple  cycle  angular  velocity  transients  were  generated  using  a  signal  or  waveform 
generator  as  a  control  input  to  the  centrifuge.  Using  the  signal  generator  as  the  centrifuge 
control  input  allowed  the  peak  to  peak  values  and  frequency  of  the  periodic 
velocity  transient  to  be  specified  for  experimental  purposes.  For  all  experiments,  the 
angular  velocity  was  increased  linearly  to  a  peak  value  at  the  midpoint  of  the  cycle  and 
subsequently  decreased  as  discussed  in  chapters  T-TTT. 

A  tri-axis  accelerometer  (Columbia  Research  Laboratories,  Inc.)  was  used  to 
monitor  the  time  variant  acceleration  components,  in  a  cartesian  reference  frame  affixed 
to  the  centrifuge  table,  in  response  to  the  angular  velocity  transient.  The  accelerometer 
output  reflects  the  acceleration  con^onents  as  referenced  to  the  mounting  location  of  the 
accelerometer  on  the  centrifuge.  A  transformation  of  the  accelerometer  output  from  the 
accelerometer  mounting  location  to  the  capillary  pivot  point  was  performed  using  the 

relations  (r'  =  =  0.8072m) 
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(i) 

pivot  ^  ^y,accel^ 

(ii) 
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(iii) 
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dt 

d(i> 

dt  >  -^r.  accei  ~  0.9819  777. 

^T,  pivot  ^ 

«  accei 

•‘•T,  accei 


Capillary  Test  Cell 

A  sealed  test  cell  was  fabricated  of  either  acrylic  or  delrin  to  contain  the  working 
fluid  and  capillary  tube  as  shown  schematically  in  Figs.  1.3,  2.4,  and  3.2  or  as  a  detailed 
drawing  in  Fig.  B.l.  The  photographs  in  Figs.  B.2  and  B.3  show  the  test  cell  as 
fabricated  of  delrin.  The  capillary  tube  was  held  in  place  using  spring  loaded  positioning 
clanps.  A  transparent  window  allowed  viewing  of  the  capillary  tube  during  experimental 
testing.  The  sealed  test  cell  contairung  a  glass  capillary  tube  and  reservoir  was  positioned 
on  the  centrifuge  as  discussed  in  chapters  I-DI.  The  test  cell  was  mounted  to  a  motorized 
optics  rotation  stage  (Newport),  as  shown  in  Fig.  B.4,  such  that  the  capillary  tube  was 
allowed  to  pivot  about  its  center  of  mass.  This  pivot  point  was  displaced  vertically  from 
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the  centrifuge  surface  at  a  fixed  radial  location,  r'o,  as  discussed  in  chapters  T-TTT.  Angles 

of  inclination  of  the  capillary  tube  were  determined  from  the  output  of  a  calibrated 
angular  displacement  transducer  (Trans-Tek,  Inc.)  accurate  to  within  ±  0.2  degrees. 

Figure  B.5  shows  the  test  cell  mounted  on  the  centrifuge  along  with  the  video 
camera  to  monitor  the  dynamic  response  of  die  meniscus.  Experiments  were  conducted 
by  fixing  various  size  diameters  of  capillary  tubes  at  a  specified  inclination  angle  while 
subjecting  it  to  the  transient  acceleration  field.  Increasing  the  inclination  angle  elevates 
the  transverse  acceleration  component  referenced  to  the  capiUaiy  tube.  This  allowed  the 
dynamic  response  of  the  meniscus  subjected  to  an  increasing  transverse  acceleration 
component  to  be  observed  experimentally.  The  dynamic  response  of  the  meniscus  to  the 
transient  acceleration  field  was  observed  using  an  8-mm  format  CCD  camera  (Sony,  30 
firames/s)  mounted  adjacent  to  the  test  cell.  The  meniscus  height  relative  to  the  reservoir 
meniscus,  t]  ,  was  measured  as  a  function  of  time  to  within  ±1.0  mm  using  the  magnified 
video  image  with  a  calibrated  reference  contained  within  the  test  cell.  An  equivalent 
contact  angle,  0 ,  was  calculated  from  the  experimentally  determined  static  wicking  height, 
hg,  for  conditions  when  the  transverse  acceleration  component  was  zero 
(hgpg  (2o/a)cose). 

Water  (reverse  osmosis)  and  ethyl  alcohol  were  used  as  the  test  fluids  to 
demonstrate  the  importance  of  the  equilibrium  contact  angle  on  the  advance  and  recession 
of  the  meniscus  due  to  a  transient  acceleration  field.  Properties  of  these  fluids  were 
assumed  to  be  that  of  the  bulk  fluid  corresponding  to  a  experimental  temperature  of  30°C. 

The  test  ceU  and  capillary  tube  were  cleaned  using  potassium  hydroxide  and  rinsed 
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with  distilled  water.  The  test  cell  was  filled  and  subsequently  sealed  resulting  in  a  fill 


that  had  a  combination  of  air,  fluid  vapor,  and  liquid. 


Experimental  Instrumentation 

Experimental  instrumentation  including,  calibration  instrumentation,  used  for  this 
investigation  are  as  follows: 


Calibration  Devices  (secondary  standards) 

(1)  CLINOMETER:  High  Resolution  (P/N  02543-01), 

Anglestar  Protractor(Schaevitz), 

Lucas  Sensing  Systems,  Inc. 

LINEAR  RANGE:  ±  19.99  deg. 

RESOLUTION:  0.01  deg. 

CALIBRATION:  USAF  Calibration  Laboratory 

Traceable  to  NIST 

(2)  TEMPERATURE 

PROBE:  Hart  Scientific  5690 

SENSOR 

ELEMENT:  High  Precision  4-Wite  Platinum  Sensor 

RESOLUTION:  O.OOrC 

TEMPERATURE 

RANGE:  -183  -  500“C 

CALIBRATION:  Factory  Calibration 

5  mK  Uncertainty,  Traceable  to  NIST 


Instrumentation  Devices 

(1)  ACCELEROMETER:  High  Precision  (SA-307-HPTX), 

Columbia  Research  Laboratories,  Inc. 
RANGE:  ±1.5g 

^3.,  ±1.0g 

OUTPUT:  ±7.5  Volts  into  lOOK  Load 

RESOLUTION:  0.001%  F.R. 

NON-LINEARITY:  ±0.03%  F.R. 

NATURAL 
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FE^QUENCY;  75  to  300  Hz  Dependent  Upon  Range 

CALIBRATION;  Factory  Calibration,  Verification  Using  Digital  Tachometer 

(2)  ANGULAR 
DISPLACEMENT 

TRANSDUCER:  Model  0605-S41020202, 

Trans-Tek,  Inc. 

RANGE;  0  -  90  deg. 

SENSITIVITY:  0.0074  VDC  /  deg.  /  Volt  Input 

CALIBRATION:  Factory  Calibration, 

Calibration  With  Assembled  Test  Cell 
Performed  Using  Clinometer 

(3)  THERMOCOUPLES:  Type-T 

SIGNAL 

CONDITIONER 

OUTPUT:  ±5  Volts 

RANGE:  0  -  199°C 

CALIBRATION:  Calibration  With  High  Precision 

Temperature  Probe  Using  a  Constant 
Temperature  Oil  Bath 


Calibration  Procedures 

Calibrations  were  performed  on  the  angular  displacement  transducer  and 
thermocouples  using  the  high  resolution  clinometer  and  the  high  precision  temperature 
probe  respectively.  Measurements  of  the  capillary  tube  diameters  were  verified  using  pin 
gauges  to  within  0.001  inch. 

Angular  displacement  transducer  use  and  calibration.  The  angular  displacement 
transducer  was  used  to  position  the  cjqpillary  tube  at  a  specified  inclination  angle  with 
respect  to  the  centrifuge  surface.  Since  the  angular  displacement  transducer  was 
calibrated  prior  to  receipt,  it  was  only  necessary  to  zero  the  transducer  with  respect  to  the 
capillaiy  tube  and  spot  verify  the  factory  calibration.  Figure  B.6  shows  the  results  of 
zeroing  the  angular  position  transducer  and  subsequent  calibration  fitted  to  a  linear 


regression. 


Thermocouple  use  and  calibration.  Themiocouples  are  a  &st  order  response  temperature 
measurement  device  for  which  dissimilar  metals  generate  a  voltage  (Seebeck  voltage). 
To  measure  temperature  accurately,  a  reference  junction  must  be  used  in  conjunction  with 
the  thermocouple.  For  most  applications,  reference  junctions  are  incorporated  within  the 
electronics  used  for  linearization  and  amplification  of  the  thermocouple  output.  Using 
thermocouples  on  rotating  devices,  such  as  with  the  centrifuge,  requires  either  the  routing 
of  the  thermocouple  wires  through  a  rotating  slip  ring  or  conditioning  the  thermocouple 
signal  prior  to  the  slip  ring.  The  former  can  result  in  significant  error  in  temperature 
measurements  due  to  the  creation  of  additional  dissimilar  metal  junctions  with  the  slip 
ring  materials.  For  these  experiments,  the  thermocouple  signals  were  electronically 
conditioned  on  the  centrifuge  resulting  in  a  ±  5  volt  ouq)ut  signal. 

Exposed  junction  copper-constantan  type  thermocouples  were  attached  to  the 
exterior  surface  of  the  c^illary  tube  with  a  graphite  slurry  to  maintain  a  thermal  bond 
with  an  outside  coating  of  high  temperature  ceramic  to  maintain  strength  integrity.  The 
capillary  tube  was  mounted  into  the  test  cell  and  the  entire  test  cell  immersed  in  a 
constant  temperature  oil  bath  for  calibration.  The  data  acquisition  system  (including 
signal  conditioner,  sl^  ring,  A/D  converter,  etc.)  and  thermocouples  were  calibrated  over 
a  temperature  range  of  25  -  70°C  using  the  high  precision  temperature  probe  as  a 
reference.  Calibration  data  for  each  thermocouple  were  obtained  and  fit  to  a  linear 
regression  for  experimental  use.  Fig.  B.7. 
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Uncertainty  Analysis 


Each  of  the  above  mentioned  instrumentation  devices  and  experimental 
measurements  have  some  error  associated  with  their  use.  This  error  must  be  considered 
either  during  the  calibration  procedure  or  when  using  instrumentation  to  determine  the 
overall  calibration  and  operational  error  due  to  the  instrumentation.  Instrument  precision, 
accuracy,  and  resolution  must  all  be  considered  when  evaluating  instrument  error. 

If  a  function,  O ,  is  considered  in  which  =  ®(a,  p,  Y),the  absolute  error 
for  $  can  be  determined^  by, 


where  A  (  )  is  the  error  associated  with  the  measurement  of  the  appropriate  independent 

variable.  When  the  A  (  ) ’s  are  considered  not  as  absolute  limits  of  error,  but  rather  as 

statistical  bounds  such  as  ±  3  s  limits,  probable  errors,  or  uncertainties,  the  proper  method 
of  combining  the  errors  is  to  use  the  root-sum  square  so  that  the  absolute  error  becomes 


Error  in  inclination  ande  measurement.  Sources  of  error  in  measurement  of  the 
inclination  angle  include  error  due  to  centrifuge  table  runout  and  voltage  measurement 
resolution.  The  error  associated  with  the  centrifuge  table  runout  was  determined  to  be 
±0.05®  and  was  added  to  the  absolute  error  associated  with  uncertainties  due  to 
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measurement  resolution.  From  the  calibration  linear  regression  and  assuming  a  sensitivity 
of  9  deg./volt  and  the  voltage  measurement  resolution  of  ±0.01  volts,  the  absolute  error 
in  inclination  angle  (using  either  error  equation)  due  to  measurement  resolution  wasE^ 
=  ±0.09°.  When  this  error  is  added  to  the  variation  of  inclination  angle  due  to  the 
centrifuge  table  runout,  the  total  absolute  error  becomes  <  ±  0.2°. 

Error  in  temperature  measurement.  Thermocouple  outputs  are  signal  conditioned  using 
a  signal  conditioner  with  a  sensitivity  of  40°CA^olt.  The  thermocouple  signal  was 
digitized  using  a  16  bit  A/D  converter  over  a  ±5  volt  signal.  This  translates  to  a  digital 
voltage  resolution  of 

±5  voi  t  _  10  voi  t 
216  2^® 

=  1, 53  xlO-^  Wit 

Mult^lying  this  digital  voltage  resolution  times  the  sensitivity  will  result  in  a 
measurement  resolution  of  6.12xlO'^°C.  The  absolute  error  in  temperature 
measurement  is  result  of  the  variability  observed  during  calibration  due  to  variations  in 
ice  point  drift  and  thermocouple  thermal  equilibrium  with  the  high  precision  temperature 
probe  during  the  temperature  bath  transients.  The  predicted  temperature  using  the  linear 
regression  fit  from  the  calibration  data  showed  a  deviation  of  less  than  TC  from  the 
actual  tenqjerature  indicated  by  the  high  precision  temperature  probe. 

Error  in  acceleration  measurement.  Accelerometer  outputs  for  the  three  axis  have  the 
following  sensitivities 
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by.  Ky 


by. 


Ky 


0.1333 
1.146  - 
0.1644 


9 

volt 

9 

vol  t 
9 

vol  t  ' 


Each  accelerometer  signal  in  digitized  using  a  32  bit  A/D  converter  over  a  ±  10  volt 

signal.  This  translates  to  a  digital  voltage  resolution  of 

±10  vol  t  _  2  0  vol  t 
216  2^2 

=  4.66  xl0‘®  volt. 

The  absolute  error  due  to  the  measurement  voltage  resolution  is  determined  by 
multiplying  the  sensitivity  for  each  axis  by  the  digital  voltage  resolution.  It  follows  that 
the  measurement  resolution  for  each  axis  becomes 

by.  Ey  =  6 .2x10-^° g, 
by^,  =  5.3x10-5^, 
by.  Ey  =  l  .7  XlO'^^g. 

The  factory  calibration  showed  an  absolute  error  in  sensitivity  would  result  in  a  less  than 
±0.013  g  error  in  accuracy.  Additional  errors  due  to  case  alignment  and  radial  position 
on  the  centrifuge  also  contribute  to  the  accelerometer  error.  Gear  train  backlash  and 
centrifuge  runout  contribute  to  transient  accelerations  during  the  experiments.  Combining 
all  of  these  potential  sources  of  error  showed  a  total  absolute  error  and  acceleration 
variation  during  the  experiments  to  be  on  the  order  of  ±0.1  g. 
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Hrrnr  in  capillary  tube  diameter.  Three  sizes  of  capillary  tube  diameters  were  used  during 
the  course  of  this  investigation.  Pin  gauges  were  used  to  determine  the  capillary  tube 
diameter  for  each  capillary  tube; 

(i) 

2a  =  1 . 892  ±0 . 038/nm, 

(ii) 

2a  =  1 . 016  ±0 . 025/nm, 

(iii) 

2a  =  0 .483  ±0 . 025/nn7. 

Error  in  video  and  time  measurement.  The  non-dimensional  meniscus  position  and  time 
were  extracted  from  video  images  recorded  on  an  8-mm,  30  frames/sec  formatted  t£5)e. 
Meniscus  height  measurements  were  made  at  specific  time  intervals  during  the  course  of 
the  acceleration  transient.  Time  intervals  for  meniscus  height  measurements  were 
specified  to  be  less  than  2%  of  the  cycle  period  for  accurate  resolution  of  the  transient 
variation  in  meniscus  height.  Meniscus  height  measurements  were  referenced  to  a  tangent 
line  to  the  bottom  of  the  reservoir  meniscus  as  shown  in  Figure  B.8.  The  absolute  error 
in  the  non-dimensional  meniscus  position  is  defined  by 


At)  ,  Aha  ~  0 . 3  c/n 
ha  ~  20  cm 
consider  two  cases; 

(i) 


n-ha, 


\2 

/ 


/ 

\ 


2 


=  0.021 


(ii) 


Tl-0, 


N 


Ati 


0.015 


Experimental  Procedures 

Experimental  parameters  were  specified  such  that  the  capillary  and  Bond  numbers 
were  of  an  order  of  magnitude  sufficient  to  investigate  viscous,  inertial,  and  transverse 
acceleration  effects  on  the  dynamic  response  of  the  meniscus  subject  to  a  transient 
acceleration  field.  These  e^rimental  parameters  consisted  of  the  radial  location  of  the 

capillary  tube  pivot  point,  Tq  =  0 . 8072 in,  static  contact  angle,  0  fluid  properties  of 

water  and  ethyl  alcohol,  p,  p,  a,  inclination  angle,  (j),  capillary  tube  radius,  a,  and 
acceleration  transient  frequency,  f.  Table  B.l  summarizes  the  experimental  tests 
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conducted  during  the  course  of  this  investigation  and  shows  the  dimensional  experimental 
parameters  and  appropriate  groupings  of  dimensionless  numbers. 

The  static  contact  angle  was  inferred  from  experimental  data  using  the  relation 

hQpsrsin<j)  -<  -^cos0^ 

3. 


Experimental  data  were  obtained  when  the  transverse  acceleration  component  was  at  a 
minimum  as  well  as  the  static  case  when  only  the  gravitational  acceleration  was  acting 
on  the  fluid  column  inclined  at  angle,  <j>.  All  meniscus  height  measurements  used  to 
evaluate  the  static  contact  angle,  0  g,  were  referenced  to  the  center  support  block  of  the 
test  cell  in  order  to  obtain  the  dimensional  meniscus  height.  Prior  to  assembling  the  test 
cell,  the  center  support  block  was  measured  using  a  micrometer  to  within  ±  0.03  mm. 
The  apparent  static  wicking  height  and  the  acceleration  term  were  used  to  calculate  the 
contact  angle.  Error  in  the  contact  angle  was  determined  statistically  between  multiple 
samples.  Table  B.2  summarizes  the  experimental  tests  and  the  appropriate  dimension  of 
the  center  support  block  and  static  contact  angle. 

As  discussed  in  Appendix  A,  analytical  solutions  are  insensitive  to  changes  in  the 
static  contact  angle  but  are  sensitive  to  changes  in  the  dynamic  contact  angle  differing 
from  the  static  contact  angle. 

Fluid  properties  for  water  and  ethyl  alcohol  used  for  this  investigation  are  shown 
in  Table  B.3.  For  experimental  tests  with  the  addition  of  heat  to  the  capillary  tube,  the 
fluid  properties  were  taken  to  be  that  of  the  bulk  reservoir  also  at  30°C. 
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Table  B.2 

Summary  of  experimental  reference  dimensions  for  analysis  of 
video  information  and  experimentally  inferred  static  contact  angle. 


exp#’s 

fluid 

reference,  mm 

deg- 

heat^  deg. 

1-5 

HP 

6.88  ±  0.01 

22.2  ±  6.6 

6-10 

QHP 

6.248  ±  0.025 

16.7  ±  4.4 

11-21 

6.713  ±  0.015 

22-24 

6.573  ±  0.025 

35.6  ±  2.0 

Table  B.3 

Summary  of  fluid  properties  used  for  this  investigation. 


fluid 

p,  kgAn^ 

H,  kg/hi-s 

0,  N/m 

HP,  30"C 

997.0 

855e-6 

0.0717 

QHgO,  30“C 

781.0 

0.001 

0.02189 

Experimental  tests  were  conducted  with  and  without  the  addition  of  heat  to  the 
capillary  tube.  In  general,  even  for  tests  in  which  heat  was  not  added  to  the  capillary 
tube,  the  acceleration  transient  was  not  applied  until  thermal  equilibrium  was  obtained. 
A  simple  string  and  weight  mounted  on  the  test  cell  allowed  for  visualization  of  the 
acceleration  vector  direction  and  peak  magnitude  in  order  to  correlate  the  video  time  to 
the  acceleration  cycle  generated  by  the  waveform  generator. 
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Fig.  B.5  Photograph  of  experimental  apparatus  mounted  on  centrifuge. 
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79a 


tcl  =  0.99402*(40°C/volt)*(volt)  -  0.39294 


181 


tc2  =  1 .0083*(40°C/volt)*(volt)  -  2.1  602 


Fig.  B.7(contd)  b)  tc2, 
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tc3  =  0.99255*(40°C/volt)*(volt)  -  0.771 1  5 
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tc4  =  0.99249*(40°C/volt)*(volt)  -  0.57348 
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tc5  =  0.99208*(40°C/volt)*(volt)  -  0.92477 


185 


tc6  =  0.991  6*(40°C/volt)*(volt)  -  0.1  6203 


Fig.  B.7(contd)  f)  tc6. 
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temperature, 


tc7  =  0.99583*(40Vvolt)*(volt)  -  0.9948 
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Fig.  B.8  Meniscus  height  measurement  reference  specification. 
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Table  B.l 

Summary  of  experimental  tests  conducted  for  this  investigation. 
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<l> 


20°  0 


30 

o 

40 

o 

50 

o 

60 

o 

25 

o 

30 

o 

40 

o 

50 

o 

60 

o 

30 

o 

40 

o 

50 

o 

60° 


25°  0 


30  ®  0 


30®  1.0 

3.5 


0.0015  0.946 


0.100  0.50 

0.050 

0.025 


0.100  0.25 

0.050 

0.025 

0.010 


0.0015  0.50 


Bo 

Ca 

2.05 

1.50x10-^ 

1.03x10-® 

7.99  xlO-’ 

6 .70x10-’ 

5.93x10-’ 

1.47 

3 . 55x10-® 

3.00x10-® 

2.33  xlO-® 

1.96x10-® 

1.73x10-® 

0.3675 

6.00x10-® 

4.67  xlO-® 

3.92x10-® 

3.47  XlO-® 

1.47 

2.37  XlO-* 

1.18x10-* 

5.92x10-® 

0.3675 

4.00x10-* 

2.00x10"* 

1.00x10"* 

4.00x10"® 

1.47 

3.00x10"® 
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exp# 

Re 

1 

0.278 

2 

0.191 

3 

0.148 

4 

0.124 

5 

0 . 110 

6 

6 .07x10-2 

7 

5.13x10-2 

8 

3 .98x10-2 

9 

3 .35x10-2 

10 

2.96x10-2 

11 

5.13x10-2 

12 

3 .99x10-2 

13 

3 .35x10-2 

14 

2.96x10-2 

15 

4.05 

16 

2.02 

17 

1.01 

00 

3.42 

19 

1.71 

20 

0 . 855 

21 

0.342 

22 

5.13x10-2 

23 

24 

2.25x10-2 


3 .30x10-2 


4.24x10-2 


5.05x10-2 


5.71x10-2 


1.93x10-2 


2.28x10-2 


2.94x10-2 


3 . 50  xlO-2 


3.96x10-2 


5.71x10-2 


7 .34x10-2 


8.75x10-2 


9 . 89  xlO-2 


1.93x10-2 


5.71x10-2 


2 .28x10-2 


a  ReCa 


ho  Bo 


4 . 58x10-® 


3 .15x10-® 


2 .44x10-® 


2.05x10-® 


2 .83x10-® 


2.39x10-® 


1.85x10-® 


1.56x10-® 


1.38x10-® 


4 .78x10-® 


3 .72x10-® 


3 .12x10-® 


2.76  xlO-® 


1.26x10-2 


3 .14x10-® 


7 .86  xlO-’ 


2.13x10-2 


5.31x10-® 


1.33x10-® 


2.13x10-’ 


2.39x10-® 


5.02x10-’ 


3 .90x10-’ 


3 .27x10-’ 


2 .42x10-® 


2.04x10-® 


1.59x10-® 


1.33x10-® 


1.18x10-® 


1.63x10-2 


1.27x10-2 


1.07x10-2 


9 .43x10-® 


1.61x10-* 


8.05x10-2 


4 .03x10-2 


1.09x10-2 


5 .44x10-* 


2 .72x10"* 


1.09x10-* 


2 .04x10-® 


192 


Table  B.4 

An  example  of  the  raw  data  extracted  from  video  information. 
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